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I.  INTRODUCTION 

The  computer  code  ATLES  and  its  associated  subprograms  and 
data  bases  are  designed  to  perform  band  model  calculations  of  the 
Infrared  radiance  signatures  of  missile  and  aircraft  exhaust  plumes  at 
arbitrary  altitude  as  viewed  through  arbitrary  atmospheric  slant  paths. 
The  code  Is  significantly  Improved  over  previous  codes  through  the  use 
of  newly-developed  band  radiation  formulations  for  hlghly-nonunlform 
optical  paths  and  the  use  of  temperature-consistent  band  model  para- 
meters for  H^O  and  CO^- 

The  band  radiation  model  Is  formulated  within  the  statistical  band 
model  for  Eorentz  line  shape  with  approximate  corrections  made  for 
Doppler  broadening  effects.  Thus,  the  model  can  be  applied  to  plume 
problems  Involving  arbitrary  target  and  sensor  altitudes.  An  optical 
line  of  sight  from  a sensor,  through  the  absorbing  atmosphere,  and 
through  a chord  of  an  exhaust  plume  Is  handled  as  a coupled  optical  path 
so  that  lln  -■  correlation  between  the  plume-emlsslon  and  atmospheric- 
absorption  spectra  Is  correctly  treated.  Line  correlation  Is  Important 
whenever  the  plume  and  atmosphere  contain  molecv  lar  species  In  com- 
mon (e.g.,  H^O  or  CO^)-  Optical  path  nonuniformities  caused  by  this 
coupling  approach  are  generally  too  large  to  be  handled  accurately  by 
the  traditional  Curtis -Godson  approximation.  For  such  coupled  paths, 
more  sophisticated  approximations  have  been  developed  and  are  employed 
In  ATLES.  Complete  descriptions  of  the  radiation  models  are  given  in 
Section  II. 

In  its  most  general  computational  mode,  ATLES  computes  the 
spectral  and  spatial  distributions  of  radiance  (both  unattenuated  and 
attenuated  by  the  Intervening  atmosphere)  over  the  projected  area  of  a 
plume  for  arbitrary  viewing  aspect  angle.  Calculation  options  allow 
integration  of  the  radiance  distributions  over  spectral  bandpasses  to 


"ATLES  is  an  acronym  for  Atmospheric  Transmittance  to  Line  Emission 
Sources  and  is  meant  to  be  pronounced  as  atlas. 
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get  inband  radiance  values  and/or  spatial  integrations  over  the  pro- 
jected area  of  a plume  to  get  axial-station  variations  or  whole  plume 
averages.  The  ratio  of  the  attenuated  to  the  unattenuated  whole  plume 
radiance  spectra  defines  an  "effective”  atmospheric  transmittance  which, 
because  of  line  correlation  effects,  will  generally  not  be  the  same  as 
that  obtained  by  traditional  atmospheric  transniittanc e codes.  The 
latter  codes  predict  attenuation  factors  for  continuum  radiation  only, 
or  for  emission  and  absorption  by  dissimilar  species.  Although  the 
primary  emphasis  of  the  calculatlonal  routines  of  ATLES  is  directed 
toward  the  plume  - atmos  phe  re  problem,  subsidiary  calculation  modes 
can  be  selected  to  simulate  laboratory  cell  and  flame  radiation  meas- 
urements. 

Primary  inputs  to  the  code  are  the  distribution  of  pressure, 
temperature,  and  species  concentration  data  (pTc  data)  within  the 
exhaust  plume,  pTc  altitude  profiles  for  the  atmosphere,  and  band 
model  parameters  for  the  molecular  species  of  the  plume  and  atmos- 
phere. The  specification  of  pTc  data  for  the  plume  is  strictly  a user's 
responsibility.  No  data  base  of  plume  models  is  maintained  for  use  in 
ATLES.  Band  model  parameters  and  model  atmospheres  may  be  user 
specified  also;  but  data  bases  of  parameters  and  atmospheres  have 
been  constructed  specifically  for  use  with  ATLES. 

The  model  atmosphere  data  base  (Table  1)  contains  six  familiar 

(2) 

standard  atmospheres. 

The  band  model  parameter  data  base  contains  parameters  for  H^O, 
CO^,  CO,  HE,  HCl  and  H^O.  Spectral  regions  and  resolutions  for  these 
various  parameter  sets  are  given  in  Table  2.  Parameter  sets  whose 
identification  name  is  prefixed  by  NASA  were  constructed  from  the  tabu- 

(3) 

lation  of  Ludwig,  et  al.  The  prefix  LINAVE  denotes  parameter  sets 

(4) 

derived  from  the  AFGL  atmospheric  absorption  line  data  compilation. 
Treatment  of  the  plume  and  atmosphere  in  a coupled  manner  reCjuires 
that  the  band  model  parameters  used  in  the  radiation  model  be  consistent 
in  the  temperature  range  from  atmospheric  values  (^ZOO^’K)  to  gas 
combustion  {'-SOOO^K)  values.  For  H^O  and  CO^.  the  NASA  parameters 


Section 

Number 


Identification 

Name 


MCCTROPICL 


MCCMIDLATS 


MCCMIDLATW 


MCCSUBARCS 

MCCSUBARCW 


Atmospheric  Model 


Tropical 

Midlatitude  Summer 
Midlatitude  Winter 
Subarctic  Summer 


Subarctic  Winter 


6 


US1962 


U.  S.  1962  Standard 


Band  Model  Parameter  File 
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have  been  shown  to  be  inadequate  for  low-temperature  application  and 
the  LINAVE  parameters  inadequate  for  high-temperature  application. 

The  latter  inadequacy  results  from  the  lack  of  hot  line  data  in  the  AFGL 
compilation.  For  the  reverse  application,  the  two  parameter  sets  are 
quite  adequate.  The  prefix  COMB  denotes  parameter  sets  that  have  been 
constructed^^’  by  combing  the  low-temperature  LINAVE  and  high- 
temperature  NASA  parameters  in  order  to  get  temperature-consistent 
parameter  sets.  Nominal  usage  considers  H^O,  CO^,  CO,  HF,  or  HCl 
as  possible  plume  radiation  species  and  H2O,  C02»  CO,  or  N^O  as 
possible  atmospheric  absorption  species. 

This  report  is  both  a description  of  and  user's  guide  for  the  plume 
radiation  code  ATLES.  The  band  radiation  formulation  is  described  in 
Section  II,  A brief  discussion  of  the  calculation  procedures  of  ATLES 
is  given  in  Section  III.  A detailed  description  of  input  data  preparation 
is  contained  in  Section  IV.  Familiarity  with  the  contents  of  Sections  II 
and  III  on  the  description  of  the  code  is  not  crucial  to  data  preparation 
or  use  of  the  code.  Example  applications  of  the  code  are  given  in  Sec- 
tion V,  and  information  on  distribution  of  the  code  to  interested  users 
is  given  in  Section  VI. 
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II.  RADIATION  MODEL 


A.  Introduction 


Radiation  calculations  in  ATLES  are  performed  entirely  within  the 
statistical  band  model  formulation  with  an  inverse  or  exponential-tail ed 
inverse  line  strength  distribution.  Derivations  and  discussions  of  the 
basic  models  and  approximations  that  account  for  nonuniformities  along 
optical  paths  and  Doppler  broadening  effects  are  given  in  Refs.  7 through 
10.  Only  the  equations  defining  the  models  and  methods  of  calculatioji  are 
given  here.  In  general,  the  optical  path  considered  is  a highly-nonuniform 
path  consisting  of  a cool  absorption  path  coupled  to  a hot  emission  path 
with  common  molecular  species  in  each  path  segment.  The  terms  "absorp- 
tion path"  and  "emission"  patH  are  identification  names  only,  since  any 
optical  path  both  absorbs  and  emits  radiation.  The  names  identify  the  1 
dominant  radiation  phenomenon  associated  with  each  path  in  the  context 
of  the  present  problem.  The  absorption  path  may  represent  either  an 
atmospheric  slant  path  or  laboratory  absorption  cell,  while  the  emission 
path  may  represent  a line  of  sight  through  an  exhaust  plume  or  a labora- 
tory flame  or  hot-gas  cell.  A schematic  diagram  of  a coupled  optical 
path  is  shown  in  Figure  1.  This  diagram  is  meant  to  convey  the  impres- 
sion that  the  optical  path  is  nonuniform  both  through  gradients  in  the 
absorption  and  emission  path  segments,  and  through  the  fact  that  the 
segments  are  treated  in  a coupled  manner.  Treatment  in  a coupled 
manner  is  mandatory  when  active  species  are  common  in  both  the 
absorption  and  emission  path  segments.  The  high  degree  of  line  position 
correlation  between  the  monochromatic  emission  and  absorption  spectra 
precludes  decoupling  the  total  optical  path  into  its  component  segments 
and  computing  the  radiance  spectra  at  the  sensor  position  as  the  product 
of  the  emission  path  radiance  and  absorption  path  transmittance.  Depend- 
ing on  the  degree  of  nonuniformity  within  the  emission  path  segment,  such 
a procedure  could  yield  a senMor  radiance  that  is  larger  or  smaller  than 
the  true  radiance,  and  only  fortuitously  the  same. 
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The  following  radiation  band  models  are  applicable  to  any  non- 

uniform  optical  path  along  which  the  pressure  p(s),  temperature  T(s), 

and  concentration  of  active  species  c(s)  is  known  as  a function  of  path 

position  s.  The  total  length  is  s^  and  the  sensor  is  located  at  s = 0. 

It  is  assumed  that  no  radiation  source  exists  beyond  s = s . The  band 

m 

model  parameters  k(s),  6(s),  and  y(s)  are  presumed  known  along  the 
path  through  p(s),  T(s),  and  c(s). 

B.  Radiance  Calculation 

The  mean  value  of  radiance  in  an  interval  Av  around  v is  given  by 


- { N*(s)  ds 


where  N*(s)  is  the  Planck  radiation  function  evaluated  at  temperature 
T(s)  and  wavenumber  v.  T(s)  is  the  mean  transmittance  in  Av  for  the 
path  between  s = 0 and  general  position  s caused  by  all  active  molecular 
species.  It  is  axiomatic  within  the  statistical  model  that  the  total  trans- 
mittance due  to  several  species  is  the  product  of  the  transmittances  for 
each  species.  Thus 

T(S)  = n~:(8) 
j ^ 


dT.(s) 


can  be  used  to  compute  t(s)  and  dT(s)/ds  once  the  transmittance  and 
transmittance  derivative  for  each  species  are  known.  The  remaining 
discussion  is  limited  to  one  species,  and  the  subscript  j is  suppressed. 
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Within  the  form  of  the  statistical  model  used  here,  rhe  mean 


transmittance  is  given  by 


W(s) 

T(s)  = e ' 

where  W(s)/6  is  an  equivalent  width  function.  The  transmittance 
derivative  is 


ds 


- T(s) 


1 dW(s) 

£ ds 


For  all  line  shapes  and  inhomogeneity  approximations,  the  equivalent 
width  derivative  is  given  by 

i 

where  only  the  function  y(s)  changes  with  the  model.  The  following 
sections  deal  with  this  y-function.  When  y(s)  has  been  determined  as 
a function  of  s,  and  thus  also  the  equivalent  width  derivative  by  Eq.  (6), 
the  equivalent  width  is  determined  by  integration 


s 

W(s)  _ r 1 dWfs')  j , 
6 6 ds'  ® 

0 


or,  in  the  case  of  the  CG  approximation,  by  formula  [Eq.  (18)  or  (43)]. 
Then,  Eq.  (4)  can  be  used  to  obtain  t(s),  and  Eq.  (5)  to  get  dT(s)/ds. 
Once  (s)  and  dT(s)/ds  are  known  for  all  species,  Eqs.  (2)  and  (3)  are 
used  to  obtain  the  total  transmittance  and  derivative  and  Eq.  (1)  used 


- 


u(s)  fe  (s)  ^ c(s-)  p(s')  k(s')  3(s')  ds' 


0(s)  . 21^ 


Within  the  CG  approximation,  the  equivalent  width  is  given  in  closed 
form  by 


In  the  more  sophisticated  derivative  approximation, 


y(s)  = y[nx(s),  p(s),  q(s)] 


where 


y('H,  P.q)  = - 


,^22 

1 + p z 


1 + z 


A numerical  procedure  for  computing  this  function  is  discussed  in  Ref.  (10). 
Briefly,  the  procedure  consists  of  rewriting  y(n,  p,  q)  m an  alternate  form 


X 

y('n,  p,q)  = - y(T|,  P,  1)  + ^ y(xz^^,  p,  l)ds 


that  relates  y(n,  p,  q)  to  its  value  at  q = 1,  The  integral  term  is  evaluated 
numerically  by  Gaussian  quadrature,  and  the  function  y(r,  p,  1)  is  approx- 


imated by 


-16- 


In  the  derivative  approximation,  W(s)/6  must  be  computed  according  to 

Eq.  (7). 

The  third  model  is  a hybrid  of  a model  developed  in  Ref.  (9)  that 
takes  greater  account  of  the  differences  in  properties  between  the  emission 
and  absorption  segments  of  the  total  optical  path  (but  still  treats  the  seg- 
ments as  coupled),  and  a model  developed  in  Ref.  (10)  that  treats  all  lines 
in  Av  as  though  they  vary  with  temperature  in  an  identical  manner.  This 
hybrid  model  (2 -path  mean  line  approximation)  is  applied  only  for  Sq  s s ^ 
(see  Figure  1).  For  s < Sq,  the  CG  or  derivative  approximation  is  used. 

In  this  model 


where 


y('Hi.  Pj,  P^) 


2 2 [A  + B(p^,p2)  + B(p2,Pj) 


- ^ 


Pi  Pz 


C(P^) 

2 2 

Pi  Pz 


C(P2) 

Pi  Pz 


y<^2’  pz> 


y(ri^,  Pj) 


Pi  = 


(1  - Pj  ) (1  - 02“) 

(1  - Pj)  (1  - p^) 


B(p,  p') 


(1  - pj^  p)  (1  - 0^  p) 
(1  - p)  (p'  - p) 


C(p)  = 


(1  - Pi^  P)  (1  - p/  P) 


/F  (1  - P)^ 


1 


(1  - p) 


2,,  2 2 2 - P,  ) (I  - P- 

Pi  (1  - P2  ) + P2  (I  - Pj  ) 4 


] Pj  ?2  1 

Pi  / P2  = 1 

P2  Pi  = 1 

: P2  and  Ti  ^ = 0 
: P2  and  "^2=0 


^ (3  - P)] 


y(n.  p)  = ^ 


'1  + 2r;  p + v/ 1 + 2t| 


1 / 

Pj  = ^ b + v/  b - 4c 


?2  = i b - yJh^  - 4c 


1 + 2t]j  1 + Zr\^ 

b = ^2  + r~ 

Pi  Pz 


1 + 2(t1j,  + 112) 

2 2 

Pi  Pz 


The  X,  p and  q parameters  In  Eq.  (Z4)  are  obtained  from  Eqs.  (IZ)  through 
(17)  and  Eqs.  (ZZ)  and  (Z3)  with  the  following  modifications:  (1)  for 
subscript  zero  parameters,  the  lower  and  upper  integration  limits  are 
0 and  Sq,  respectively;  (2)  for  nonsubscripted  parameters,  the  limits  are 
Sq  and  s,  respectively. 

D.  y-Functions  for  Doppler  Line  Shape 
In  the  CG  approximation 

y(s)  = y[TTx(s),  p(s)] 


y(Ti,  p)  = (2  - p)  + (p  - 1) 
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1,1  I.  Vf. 


i 


t 


and 


H(n) 


1 + Tie 


du 


(36, 


The  functions  H(ri)  and  dH(ri)/dTi  can  be  obtained  by  interpolation  on  the 
data  of  Table  3.  For  small  ti 


H(n) 


(-T1) 


n + 1 


(n  + 1) 


3/2 


dHQl) 

dri 


E 


n = 0 


v/  n + 1 


while  for  large  ri 


(37) 


(38i 


H(p)  - 

4 

3 

(In-n)^''^ 

1+ 

(InT^)'^ 

, 7it"^/640 

+ 4 + . . . 

(In-)^ 

(39' 

* 

dH(ri)  ^ 

2 

y Inri 

1 tt^/24 

(Inl])^ 

35r'^/1920 

(40 

dri 

4 + • . . 

(Irni)^ 

! 


The  arguments  x(s)  and  p(s)  are  defined  as  for  the  Lorentz  line  shape  by 
Eqs.  (12)  through  (16),  but  with  P(s)  defined  as 


0(s) 


Y(s) 

6(s) 


where  y(s)  is  the  Doppler  line  width 


(41 


I 
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Table  3.  Curve  of  Growth  FunctLons  for  Statistical 
Doppler  Line  Band  Model 


•n 

H(ti) 

dH(r|)/dri 

0.  1 

0. 09665 

0. 93460 

0.  2 

0. 18722 

0. 87828 

0.4 

0.35325 

0. 78593 

0.  6 

0.  50289 

0. 71310 

0.  8 

0. 63940 

0. 65398 

1.0 

0.  76515 

0. 60490 

2.0 

1.28138 

0. 44564 

4.0 

2. 00567 

0. 29978 

6.  0 

2.  52794 

0. 22935 

8.  0 

2. 94136 

0. 18717 

10 

3.28568 

0.  15883 

20 

4.47964 

0. 092598 

40 

5.84661 

0. 052199 

60 

6.71916 

0. 036911 

80 

7.36866 

0. 028752 

100 

7. 88910 

0. 023642 

200 

9.59377 

0. 012759 

400 

11.4233 

0. 006813 

600 

12. 5478 

0. 004702 

800 

13. 3689 

0. 003609 

1000 


14.0187 


0. 002937 


! 


Y(s)  3.56817  X 10* 

(v  has  unit  cm  ^ , T(s)  has  unit  °K  and  m is  the  atomic  mass  number  of  the 
molecular  species.)  Note  that  8{s)  for  the  Doppler  model  is  defined  in  terms 
of  the  parameter  5(s)  appropriate  for  the  Lorentz  model.  This  is  a necessary- 
approximation  since  an  adequate  definition  for  the  mean  line  spacing  para- 
meter in  the  Doppler  statistical  model  does  not  exist. 

The  equivalent  width  within  the  CG  approximation  for  the  Doppler  line 
shape  Is 


W(s)  _ ^e^®^ 


H[ttx(s)] 


The  Doppler  equivalent  of  the  Lorentz  derivative  approximation 
[Eq.  (19)],  is  not  convenient  for  practical  application.  A related  model 
based  on  the  mean-line  approximation  discussed  above  is  used  in  this 
case.  The  y-function  is 


y(s)  = y 


iq(s)  ' " 


where 


y(n.  P)  = -I-  f 


Q jl  + r exp(- 


For  small  n,  y(r|,  p)  can  be  computed  from 


y(i^.  p)  ^ ^ 


n p;  0 \ l + n 


whereas  for  p 0.  5,  the  approximation 


y(Tl.  P) 


is  accurate  to  better  than  one  percent.  For  other  (rj,  p)  conditions,  y(T|,  p) 
can  be  obtained  for  most  practical  applications  by  interpolation  on  the 
data  of  Table  4. 


Similarly,  the  Doppler  equivalent  of  the  2-path  mean  line  approxi- 
mation [Eq.  (24)],  for  Lorentz  lines  is  not  practical.  The  following  2 -path 
approximation  is  based  on  a rational  approximation  of  the  Doppler  line 
shape  by  a parabolic  line  shape  [see  Ref.  (9)].  The  y-function  is 


y(s) 


nx(s) 

q(s)' 


^0’ 


where 


y(Tii.  '^2.’  Pi’  “ 


P2^Pl^l 


g(Tl^,  ri^,  Pj,  9^.  0,  1) 

Pj,  p^.o,  i/pj)  + g{0,x]^,  p^,  9^,  1/pj,  1) 

Pj,  P2«  1/Pj)  + g(0,  P^ , P2»  1 / P j . 1 / P2)  + gg^ P2>  ^ Pi 

g(^2»"n2^Pl’P2’^’^^  Pl^p2^^ 

Pr  Pz'®'  ^^Pz^  Pf’  Pz’  ^^Pz'  Pi^^^Pz 

g (tij , '^2  > Pj’  p2’^’  Pz  ^ ^ g(^2’^’  Pj’  Pz  ’ Pz’  Pl^  ^ ®0^Pl^  ^^Pl^Pz 


D^^lSPl 
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Table  4.  Equivalent  Width  Derivation  Function  for  Statistical  Doppler  Line  Band  Model 


w 

o 

o 


W 

o 

o 

(NJ 


W 

o 

iD 


|o 

I w 

1 

(VJ 


w 

o 

o 


I 

I w 

o 

o 


o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

K 

w 

w 

w 

w 

w 

w 

w 

w 

w 

w 

w 

W 

w 

w 

w 

w 

w 

r.-; 

w 

w 

w 

m 

m 

CO 

oc 

O 

03 

00 

03 

o 

r— 1 

03 

in 

CO 

o 

vD 

o 

00 

•o 

CO 

o 

r-i 

m 

vO 

O 

sD 

CO 

ON 

sO 

(— « 

CO 

o- 

ON 

03 

ON 

c:^ 

r-H 

00 

vO 

o 

. -» 

in 

nO 

1— < 

nD 

CO 

n£) 

03 

sO 

in 

CO 

CO 

i-s 

rO 

in 

o 

CT' 

CT' 

00 

o 

03 

00 

03 

ON 

ON 

m 

03 

00 

sD 

r— < 

o 

fN- 

in 

O' 

CO 

CO 

O' 

O 

00 

CO 

o- 

sO 

m 

in 

in 

CO 

CO 

CO 

ro 

o 

O 

o 

CD 

o 

O 

o 

o 

o 

CD 

O 

o 

CD 

o 

CD 

O 

o' 

o' 

o 

o* 

o’ 

o’ 

o 

O 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

WWW 

WWW 

w 

w 

w 

WWW 

w w 

r t 
U-i 

w 

w 

WWW 

w 

w 

sO 

r- 

O 

CO 

03 

00 

03  nO  nO 

O 

,-H 

cc 

00 

o 

O 

vC 

00 

o 

CO 

00 

o- 

in 

sO 

r- 

CO 

o 

O 

o 

m 

o 

m 

00 

CO 

o 

03 

rH 

m 

in 

o 

CO 

o 

ro 

vO 

O'  00 

m 

CO 

O CO  CO 

r-' 

CO 

1-^ 

03 

o 

03 

m 

00 

CO 

ON 

m 

03 

O 

r- 

m 

03 

O'  O O' 

O O'  00 

00 

m 

CO 

CO 

03 

03 

03 

03 

r-H 

o 

o 

o 

o 

o 

CD 

o 

o 

CD 

CD 

CD 

CD 

o 

O 

CD 

o' 

CD 

CD 

o' 

o’ 

o’ 

o’ 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

l-H 

1^ 

r-c 

w 

wwwwwwwwwwwww 

w 

w w 

w w 

w 

w w 

sO 

(J^  sO 

o 

CO 

03 

o 

o 

03 

m 

o- 

o cn 

1-^ 

03 

o 

sO 

in 

<-H 

CO 

o 

in 

o 

m 

in 

in 

On  O 

CO 

o 

m 

r- 

03 

00 

03 

o 

m 

p*H 

00 

00 

03 

in 

r- 

03 

in 

00 

O'  on 

CO 

03 

CJN 

o 

o r- 

m 

03 

CO 

m 

t-H 

v£>  sO 

CO 

o o 

03 

in 

o 

m 

o 

CO 

CO 

CO 

Cjs  <js  00 

0- 

CO 

03 

03 

ON  vsD 

CO 

fO 

03 

P—1 

o 

o 

o 

o 

O 

o 

CD 

o 

O 

CD 

CD 

o o 

O 

o 

CD 

o' 

d 

CD 

o 

CD 

CD 

o o 

w w 

o 

un  o 

^ CT' 
00 


o o 

■+•  4- 

w w 

vO 

(M 

rn  00 

(JN  On 


o o o o 
W W w w 

r--  O'  o 
^ cO  O 
CO  r--  0^  ^ 
O (G'  rf 

o <r\ 


o o o 

•f-  ■+•  + 

WWW 

CO  vO  CO 

00  to 
10  0^0 
sO  LO  00 
CO  <vl  «-i 


w w w w w w 

tn  Tj’  CO  (\)  03 

00  00  O O (\3 

o CO  h-  o fo 

o in  ^ m 

1-H  r-  CO  fsj 


03  OJ  03  03  rv3 

r I I I I 

M U W H W 

O ro  00  fO  O 

Tf'.osooor^ 

o CO  o ro 

LO  c~  o >-< 

O ^ r<^  CO 


oooooooooooooooooooooo 


OOOOOOOOOOOr-<~<,-<^ 
+ + 4'  + H“4-  + -h  + + + < * I I 

wwwwwwwwwwwwwww 

ooino-ocoo3inf\3coinooo3COc^ 
ooo-ooG'cocoinoorfcr^cNjTfT^ir^ 
cnr^Ocsi*— ‘'^co^*— isO*— •r^mvDoo 
COOONOTfC^r^rfO^OCOOcOsOcO 
O^CT'O^O^OOr^sOTfcOOJ'— 


03  03  03  fvj  03  CO  fO 
I I r I I I I 

w W U K W w w 

Tf  00  CO  CO  ■rf' 

O CO  rO  O CO  4* 

IT)  CT'  in  CO  o 

vf  _i  so  so  o m CO 

CO  m co  ■-H  ^ in  m 


odd 

o' 

o’ 

o' 

o’ 

o' 

o 

O 

CD  O 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

o o o 

o 

o 

o 

o 

o 

o 

O 

t-H  i-H 

.—1 

p— 1 

(N3 

rv3 

03 

CO 

CO 

CO 

1 

4-  4-  + 

4- 

-I- 

4- 

4“ 

4- 

4- 

4- 

i 1 

1 

1 

1 

1 

1 

1 

t 

1 

1 

wwwwwwwwwwwwwww 

w w w w w 

w w 

CO  00  CO 

o 

00 

o 

O 

in 

CO 

00  CO 

03 

ON  nD 

CO 

03 

CO 

sO 

O 

rO 

O'  o o 

nO 

03 

m 

ON  sO 

00 

On  n- 

in  nD 

sO 

CO 

in 

o 

sO 

03  vO  O 

00 

'if 

in 

ON  00 

ll'>  p-^ 

in 

CO 

o 

sO 

ON  -H 

C7N  00  sO 

CO 

(N- 

o 

in 

in 

03  CO 

CO 

fv3 

m 

o 

m 

in 

n- 

CO 

00 

QN  On  ON  <DN  CO 

r- 

vO 

<N3 

CV3 

in 

<N3 

vO 

CO 

CO 

o o o 

o 

o 

CD 

o' 

o 

CD 

CD 

d d 

d 

d 

o’ 

d 

d 

d 

d 

d 

d 

d 

cococommm-^-^rf 


OOOOOOOOOOrt. 

WWWWWWWWWWWUWWWWUWKWWfcl 
^roc-iroOOoOO'OcO-«~i)Tj.oOrrcOcOOooS)^fnf^ 
'-^cfcoO'Orf<vOrocooco'^ooococnooo\rcca'rti 
'^oir^ooooo3003'>ocn'*^in'^rnoLnr^oinr^OLn 
2^00ir)F-tTtc^c0s000O^o3F-,tN.cn^0.r0F-^r-mF-H 
G^CT^o^G^oosOLncn«— 4,-huSo4>— <ino>3>— <irio3F-4UD(\3^ 


oooooooooooooooo 


o o o o o o 


03  03  03 


fT)  co  m 'jf 


^ o o o -• 

HWWWWWUWWWWWWWUirjUWtijWHU 
oooooooooooooooSo''  — — 


OJ  03  03 

-F 


lTi 
-r  + 


03  m 


03  m 


03  in 


03  m 


o 

• ■ • • 

03  m ^ 03  in 


o o o o o 


i—t  03  in  *-H 


g(f'i . "fi,.  Pi . P,.  a,  b)  = y 


b-a  , (x-O')  , cv  + b V c^x  o^-av/ckx 

— + -i— In — y— 

2x  V cfx  n-b  V (vx  n + a v Q'x 


go(p)  = 


(p^  - 2p  - 2) 

2p^ 


O'  = 1 + Tj  + Tl^ 
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The  X,  p,  and  q arguments  of  Eq.  (48)  are  computed  with  the  same  Integra- 
tion limit  modifications  as  are  used  In  the  Lorentz  2-path  model.  3(s)  Is 
computed  from  Eqs.  (41)  and  (42). 

E.  y-Functlon  for  Voigt  Line  Shape 

The  band  model  formulation  for  Voigt  line  shape  Is  based  on  an 
approximation  devised  by  Rodgers  and  Wllllams^^^^  for  Isolated  lines  and 
uniform  optical  paths.  The  equivalent  width  and  y-functlon  are  computed 
assuming  first  a pure  Eorentz  line  shape  and  second  a pure  Doppler  (or 
parabolic)  line  shape  according  to  the  methods  of  the  preceding  paragraphs. 
The  equivalent  width  for  the  Voigt  line  Is  then  approximated  as 


W(s) 

6 


Wi^(s)  ^ Wj^(s) 

I 6 


^ ^WJs)  Wi^(s)  ^ 

- 5 5 / 


W^(s) 

^ = u(s)  k^(s) 


In  the  single-path  models  and 
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(56) 


W;^(s) 


- UQkQ  + u(s)k(s) 


in  the  2-path  models  for  s > s^. 


The  derivative  function  is  given  by 

y(s)  = A(s)  y^ls)  + B(s)  yj^(s)  + c(s) 


(57) 


where 


A(s)  = 


W(s)/6 


'"w^(s)/6^^ 


1 

\^V 


(s)/6y 


(58) 


B(s)  = 


Wj3(s)/6 

W(s)/6 


1 -I 


W^(s)/6 

W^(s)/6, 


C(s)  = 


W^(s)/6  /w^(s)/6  /W^{s)/bf 

W(s)/6  \W^{s)/bJ  ( W^(s)/6, 


(59) 
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III.  COMPUTER  CODE  AXLES 


A.  Flow  Diagram 

The  main  code  ATLES  and  its  associated  subprograms  are  coded 
in  FORTRAN  IV  compatible  with  The  Aerospace  Corporation  CDC  7600 
computer  and  SCOPE  operating  system.  A simplified  flow  diagram  for 
ATLES  is  shown  in  Figure  2.  This  section  gives  a brief  discussion  of 
each  of  the  numbered  blocks  on  the  flow  chart.  Sections  B through  E 
give  more  detailed  information  on  critical  computation  routines. 


1.  Each  run  of  the  program  requires  the  input  of  various  num- 
erical data  and  control  data  to  specify  the  operating  mode 
of  the  program.  The  preparation  of  input  data  is  detailed 
in  Section  IV.  The  read-in  of  all  data  for  a computational 
run  is  effected  by  a call-to -subroutine  INPUT.  This  sub- 
routine either  reads  data  directly  from  the  card  deck  input 
stream  or  calls  other  subroutines  that  read  data  from  the 
input  stream  or  from  attached  tape  or  disc  input  devices. 
Subroutines  that  may  be  called  by  INPUT  and  their  function 
are: 


MODATM 

BAND 

PLUME 

PDATA 

SGRID 


Read  altitude  profiles  of  a model  atmosphere. 

Read  band  model  parameters  for  a single  species. 

Read  pTc  data  for  an  exhaust  plume. 

Read  pTc  data  for  a single  line  of  sight  through 
an  absorption  or  emission  cell. 

Read  plume  projected  area  integration  coordinates. 


Data  read  directly  by  INPUT  include: 


a.  Control  data  specifying:  (1)  Whether  calculations  are 
to  be  made  on  an  isolated  absorption  path,  an  isolated 
emission  path  or  a coupled  absorption-emission  path; 
(2)  Whether  the  absorption  path  is  an  atmospheric 
slant  path  or  an  absorption  cell  path;  and  (3)  Whether 
the  radiation  source  is  a plume  or  an  emission  cell 
path. 
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Figure  2.  ATLES  Flow  Diagram  (Continued) 
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b.  Spectral  data  defining:  (1)  The  spectral  region  over 
which  calculations  will  be  made  and  the  step  size; 

(2)  The  spectral  unit;  and  (3)  Bandpasses  over  which 
computed  spectral  data  will  be  averaged. 

c.  Band  model  data  specifying  the  lineshape  and  inhomo- 
geneity modes  to  be  used. 

d.  Integration  step  size  data  for  path  integrations  over 
the  absorption  and  emission  portions  of  the  optical 
line  of  sight  and  for  area  integration  over  the  projected 
area  of  a plume.  The  latter  data  consist  of  grid 
coordinates  defining  many  lines  of  sight  through  the 
plume . 

e.  Information  specifying  those  species  to  be  included  in 
the  absorption  path  and  those  to  be  included  in  the  emis- 
sion path. 

f.  Geometry  data  defining  the  atmospheric  slant  path  and 
plume  aspect  angle. 

2.*  If  absorption  over  an  atmospheric  slant  path  is  to  be  com- 
puted, calls  to  the  subroutines  VIEW  and  APATH  are  made. 
The  former  computes  the  altitude  variation  over  the  slant 
path  as  a function  of  distance  along  the  path,  and  the  second 
interpolates  for  pTc  data  along  the  path  using  the  altitude 
variation  determined  by  VIEW  and  the  model  atmosphere 
read  in  by  MODATM.  The  number  of  points  along  the  slant 
path  and  the  geometry  of  the  slant  path  are  defined  as  part 
of  the  input  (step  1).  If  an  absorption  cell  is  being  considered, 
the  pTc  data  along  the  path  are  computed  in  subroutine  LPATH 
by  linear  interpolation  on  data  read  in  by  PDATA. 


*More  detailed  information  on  subroutines  MODATM  and  VIEW  is  given 
in  Sections  B and  C,  respectively. 
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3.  A spectral  loop  over  limits  and  increment  adined  as  part  of 
the  primary  input  is  now  initiated.  If  the  probleni  consists  ol 
computing  source  radiation  only,  a skip  to  step  8 occurs 
now. 

4.  Band  model  calculations  are  performed  for  one  species  at  a 
time.  This  step  initiates  a loop  over  all  species  specified 

by  the  primary  input  data  to  be  in  the  atmosphere  (or  absorp- 
tion cell). 

5.  Subroutine  BPARAM  is  used  to  compute  band  model  para- 
meters along  the  absorption  path  with  the  data  read  in  by 
RAND  and  the  pTc  variation  determined  in  step  Z.  Linear 
interpolation  on  both  v and  T is  employed  except  for  T out- 
side the  range  of  the  input  band  model  parameters  table, 
Semilogarithmic  extrapolation  is  used  in  this  case. 

6.  ' Subroutine  DELTRAN  is  the  core  of  the  band  modeling 

routines.  With  the  pTc  and  band  model  parameters  deter- 
mined in  steps  2 and  5,  respectively,  the  transmittance  and 
transmittance  derivative  along  the  optical  path  arc  computed 
according  to  the  models  of  Section  II.  Subroutines  that  may 
be  called  by  DEI. TRAN,  their  puroose,  and  the  method  of 
computation  are: 

FI  Compute  equivalent  width  function  for  a band  of 

Lorentz  lines  (analytic  function). 

GI  Compute  equivalent  width  function  for  a band  of 

Doppler  lines  (asymptotic  series  and  quadratic, 
semilog  interpolation). 

WMIX  Compute  equivalent  width  for  a band  of  Voigt  lines 
according  to  the  approximation  of  Rodgers  and 
Williams  (analytic  function). 


More  detailed  information  on  subroutine  DEL  IRAN  is  given  in  Section  E. 
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yCGL  Compute  equivalent  width  derivative  function  for  a 
band  of  Lorentz  lines  in  the  Curtis -Godson  approxi- 
mation (analytic  function). 

YDRL  Compute  equivalent  width  derivative  function  for  a 
band  of  Lorentz  lines  in  the  derivative  approxima- 
tion. This  subfunction  uses  the  subprogram  YLSL 
which  computes  the  derivative  function  in  the 
Lindquist-Simmons  approximation  (numerical  inte- 
gration by  Gaussian  quadrature). 

Y2MLL  Compute  equivalent  width  derivative  function  for  a 

band  of  Lorentz  lines  in  the  2-path,  mean-line  approx- 
imation (analytic  functions).  YLSL  is  used  by  this 
subfunction  also. 

YCGD  Compute  equivalent  width  derivative  function  for  a 
band  of  Doppler  lines  in  the  Curtis-Godson  approxi- 
mation (asymptotic  series  and  parabolic,  semi- 
logarithmic  interpolation). 

YMLD  Compute  equivalent  width  derivative  function  for  a 
band  of  Doppler  lines  in  the  mean-line  approxima- 
tion (series  solution,  analytic  approximation  and 
tabular  interjxilation). 

Y2MLP  Compute  equivalent  width  derivative  function  for  a 
band  of  Doppler  lines  (parabolic  approximation)  in 
the  2-path,  mean-line  approximation.  The  subfunction 
FUNC  is  called  by  Y2MLP  (analytic  functions). 

YMIX  Compute  equivalent  width  derivative  function  for  a 
band  of  Voigt  lines  using  the  approximation  of 
Rodgers  and  Williams  (analytic  functions). 

When  the  transmittance  and  derivative  for  each  species  in 
the  absorption  path  has  been  determined,  the  total  trans- 
mittance and  radiance  of  the  path  is  computed  by  subroutine 
RADNCE.  The  path  integration  is  handled  by  trapezoidal 
quadrature. 
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8.  If  the  radiation  source  is  a plume,  a two-dimensional  loop 
over  the  projected  plume  area  is  now  initiated  using  the 
coordinate  data  read  in  by  SGRID.  If  the  sour'.e  is  an  emis- 
sion cell,  the  loop  is  gone  through  just  once.  If  only  an 
absorption  path  is  being  considered,  no  loop  is  made,  and 
an  immediate  continuation  of  the  spectral  loop  takes  place. 

9.  If  the  source  is  a plume,  subroutine  SFATH  is  used  to  com- 
pute the  pTc  variation  over  the  portion  of  the  line  of  sight 
within  the  plume  (the  emission  path).  The  plume  data  read 
in  by  PLUME  is  used  for  this  purpose.  If  the  source  is  an 
emission  cell,  the  pTc  variation  is  computed  by  subroutine 
LPATH  and  data  read  in  by  PDATA. 

10-12,  Steps  10  through  12  are  similar  to  steps  4 through  7 except 
that  the  total  radiance  for  the  emission  source  is  computed. 
These  steps  compute  the  source  radiance  for  no  attenuation 
by  the  absorption  path. 

13-14.  Again,  these  steps  are  similar  to  steps  10  through  12  except 

that  now  atmospheric  attenuation  is  included.  If  only  an  emis- 
sion path  is  being  considered,  this  step  is  bypassed.  Note 
that  the  band  model  parameters  over  the  emission  path  do 
not  have  to  be  computed  in  this  routine  since  they  have  al- 
ready been  computed  in  step  10. 

15.  Continue  loop  over  all  lines  of  sight  through  plume. 

16.  When  the  radiance  of  all  lines  of  sight  through  the  plume  has 
been  determined,  spatial  integrations  over  the  plume  are 
performed  in  subroutine  PLUMAVE.  Trapezoidal  quadrature 
in  both  plume  directions  are  used  for  these  integrations. 

17.  Continue  spectral  loop. 
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18.  At  the  completion  of  step  17,  the  spatial  and  spectral  varia- 
tion of  both  the  unattenuated  and  attenuated  plume  radiance 
has  been  determined,  and  the  average  spectra  of  these  two 
radiances  over  the  projected  area  of  the  plume  has  been 
computed.  Subroutine  BNDPASS  is  now  called  to  compute 
bandpass  averages  of  the  average  radiance  spectra  and  of 
atmospheric  transmittance  spectra.  Linear  interpolation 
and  trapezoidal  quadrature  are  employed  in  this  routine. 

19.  Computed  results  are  now  listed  by  subroutine  OUTPUT  and 
(optionally)  plotted  by  subroutine  SPLOT.  The  nominal  list- 
ing consists  of  a summary  listing  of  input  data;  the  spectra 
for  (1)  Atmospheric  or  absorption  path  radiance;  (2)  Un- 
attenuated source  radiance  averaged  over  plume;  (3)  Atten- 
uated source  radiance  averaged  over  plume;  (4)  Atmospheric 
or  absorption  path  transmittance;  (5)  Effective  atmospheric 
or  absorption  path  transmittance;  and  bandpass  averages  for 
each  of  these  five  spectra.  An  optional  listing  of  the  spatial 
variation  of  bandpass  averaged  radiances  over  the  projected 
area  of  a plume  is  made  by  subroutine  MAP. 

20.  Return  to  beginning  of  program  for  a new  computational  run. 


B.  Subroutine  MODATM 

This  subroutine  is  designed  around  the  model  atmosphere  profiles 

(2) 

of  McClatchey,  et  al.  ' ' Data  defining  the  pressure,  temperature  and 
water  concentration  at  the  specific  altitudes  0,  1,  2,  3,  ...,  24,  25,  30, 
35,  ..  . , 50,  70  and  100km  are  read  in.  The  input  data  deck  structure 
is  given  in  Figure  13.  The  pressure  and  water  concentration  data  may 
be  entered  in  natural  units  of  atm  and  mole  fraction,  respectively,  or 
in  the  McClatchey,  et  al.  , units  of  mb  and  g/m  , respectively.  If  the 
latter,  conversion  to  natural  units  is  made  by 


and 


c.j  _ (mole  fraction)  = 4,  559  x 10 
H^O 


Ph  o 

- 6 2 

p(atm ) 


Data  entered  at  the  altitudes  specified  above  are  used  to  construct  an 
altitude  Interpolation  table  with  the  pTc  data  given  at  each  1-km  step 
between  0 and  100km.  Linear  Interpolation  on  T,  log(p)  and  log(Cj_j^^) 
for  20  s z £ 50  and  quadratic  Interpolation  on  these  same  quantities  for 
50  s z s 100  are  used  to  construct  the  table. 


C.  Subroutine  VIEW 

Subroutine  VIEW  accepts  data  defining  the  geometry  configuration 
of  a sensor  and  target  and  computes  the  altitude  variation  along  the 
slant  path  between  these  two  positions.  Two  geometry  cases  are  con- 
sidered. The  first  consists  of  a sensor  above  the  atmosphere  and  a 
target  within  the  atmosphere.  The  data  entered  for  this  case  are  the 
target  altitude  h,  the  target-to-sensor  zenith  angle  9,  and  an  altitude 
z^  defining  the  top  of  the  atmosphere.  The  slant  path  defined  by  these 
variables  Is  shown  In  Figure  3.  The  path  position  s = 0 Is  defined  as 
the  point  where  the  line  of  sight  from  the  sensor  to  the  target  enters 
the  atmosphere.  The  total  path  length  In  the  atmosphere  Is 


/ 2 2 

s^  = - (Rg  + h)  cos  0 + z^)  - (R^  + h)  sin  9 


where  R Is  the  radius  of  the  earth.  The  value  R = 6368km  Is  used, 
e e 

For  this  geometry,  and  for  any  value  of  0 ^ s s;  s^,  the  altitude  at  s Is 
computed  from 


/ 2 2 

7.(s)  + s - 2s  + (R^^  4 h)  cos  9 
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The  maximum  zenith  angle  allowed  for  this  case  Is 


, / R ' 

^ . -1  / e 

0 = TT  - sin  r5 — — ^ 

m \ R + h 


Nominally,  z = 100km  Is  used.  For  h = z then,  the  maximum  allowed 
•'  m m 

zenith  angle  Is  ~ 100  . 

The  second  viewing  geometry  case  consists  of  a sensor  within  the 
atmosphere  and  a target  either  within  or  outside  the  atmosphere.  The 
data  entered  for  this  case  are  the  sensor  altitude  H,  the  sensor-to-target 
zenith  angle  ?p,  the  range  from  the  sensor  to  target  r,  and  again  z^.  The 
slant  path  for  this  case  and  with  the  target  within  the  atmosphere  Is  shown 
In  Figure  4.  The  target  altitude  Is  computed  by 


h = + H)^  + r^  + 2(Rg  + h)  r cos  cp  - R^ 


The  total  path  length  In  the  atmosphere  Is  computed  from 


and  the  altitude  variation  along  the  path  for  0 s s s s^  by 


h s z 


- (R  + H)  cos  cr  + \/(R  + z - (R  + H)^  sln^  cp  h > z 

e Verne  m 


z(s)  = V ^ - 2s(R^  + H)  s cos  cp  - R^ 


If  the  zenith  angle  is  greater  than  the  critical  angle 


the  line-of-sight  from  the  sensor  intersects  the  earth.  In  this  case,  the 
range  r must  be  less  than  the  critical  length 


(R^  + H)  cos  cp 


-V 


- (R  + H)^  sin^  cp 
e e 


Neither  of  the  geometry  cases  consider  line  of  sight  curvature 
due  to  refraction. 

The  optical  path  within  the  atmosphere  is  divided  into  N intervals 
(N  is  an  input  parameter  to  the  subroutine),  and  the  altitudes  at  the 
N + 1 points  between  s = 0 and  s = s^  computed  according  to  the  pre- 
vious equations.  The  division  is  made  so  that  each  interval  contains 
approximately  the  same  mass.  This  is  accomplished  by  assuming  an 
isothermal,  exponential  pressure  atmosphere,  which  gives  the  algorithm 


s^  ^In 

2 cos  a 

-(Sj/P)  cos  O'  j 

° - N 

-(s  /p)  cos  a' 

1 ^ 

1 - e 

L J 

Given  the  path  position  s^,  s^  5*  Sj  is  the  path  position  such  that  1/N  of 
the  integrated  mass  along  the  path  lies  between  s^  and  s^.  A starting 
value  of  Sj  = 0 is  used.  The  computed  value  of  s^  becomes  the  next 
value  of  Sj^,  etc.  The  process  is  repeated  N times,  p is  the  exponen- 
tial scale  height  of  the  atmosphere.  The  value  B = 7.  32km  is  used. 

For  the  geometry  case  of  a sensor  within  the  atmosphere,  a = cp.  For 
the  exoatmospheric  sensor  case,  a = rr  - 9.  For  a within  z”  of  90° 
(nearly  horizontal  paths),  a simple  equal- interval  division  of  the  path 
is  made. 
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The  final  results  of  the  subroutine  are  the  two  arrays  s(l),  z(i), 
i = 1,  N + 1.  These  arrays,  along  with  the  1 -km  interval  atmosphere 
profile  table  computed  in  subroutine  MODATM  are  passed  to  subroutine 
APATH  where  an  interpolation  calculation  is  used  to  obtain  pTc  data  at 
the  N + 1 points  along  the  atmospheric  slant  path. 

D.  Subroutine  SPATH 

Subroutine  SPATH  computes  the  geometry  of  intersection  (Figure 
5)  between  a line  of  sight  and  a paraboloid  plume  control  surface  and 
interpolates  for  pTc  data  along  the  intersection  chord.  The  plume  geo- 
metry data  (nozzle  exit  radius),  L (plume  length),  and  (maximum 
plume  radius)  are  entered  as  input  data  through  subroutine  PLUME. 
Subroutine  PLUME  also  enters  the  pTc  data  defining  the  plume  within 
the  paraboloid  plume  control  surface  (see  Figure  14  for  preparation  of 
plume  data).  The  line  of  sight  through  the  plume  is  defined  by  three 
parameters:  z^  is  the  axial  distance  at  which  the  line  of  sight  intersects 
the  y-z  plane  (see  Figure  5);  d is  the  transverse  distance  from  the  x-z 
plane  where  the  line  of  sight  enters  and  exits  the  plume  control  surface; 
and  a is  the  aspect  angle  of  the  line  of  sight  with  the  plume  axis  (0  s o 5 
180°). 

The  values  of  R , R , L,  z , d,  and  a determine  one  of  five  basic 
o m o 

intersection  geometries  that  depend  principally  on  which  of  the  three 
surfaces  (nozzle  exit  plane,  paraboloid  surface,  end  plane)  the  line  of 
sight  enters  and  exits  the  plume.  These  cases  are  illustrated  in  Figure 
6.  The  case  where  the  line  of  sight  cannot  enter  the  plume  because  of 
missile  obscuration,  or  where  the  line  of  sight  does  not  intersect  the 
plume  are  handled  as  a sixth  case  by  setting  an  obscuration  flag.  The 
cases  for  o = 0 and  0 = 7:  are  also  handled  separately  from  the  five  basic 
intersection  geometries. 

The  intersection  cases  are  determined  by  comparisons  with  the 
four  test  quantities 
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Figure  6.  General  Intersection  Cases.  NOTES: 
(1)  Cases  2 and  4 handle  the  case 
where,  for  a ~ir/2,  the  LOS  exits 
from  the  control  surface  on  the  same 
side  that  it  enters.  (2)  Cases  3 and  5 
can  occur  only  for  a>  tt/2  and 
d < Rq.  If  d > Rq,  the  LOS  continue 
to  the  control  surface  and  are  handled 
by  cases  2 and  4,  respectively. 
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The  significance  of  these  quantities  can  be  deduced  from  Figure  7.  If 
z < z . or  z > z,  (a  < n/2)  or  z > z.(a  > n/Z),  the  LOS  misses  the 
plume.  Likewise,  if  d > the  LOS  misses  the  plume.  The  test 

parameter  z^  is  relevant  only  if  d < and  is  used  to  test  for  obscura- 
tion (O'  < n/2)  or  for  cases  where  the  LOS  ends  on  the  nozzle-plane.  The 
LOS  is  obscured  if;  (1)  a < n/2  and  d < R and  - z_  < z < z_;  or  (2)  If 
d < R^  when  0=0.  The  test  parameters  z^  and  z^  are  used  to  differen- 
tiate between  LOS  that  exit  (a  < n/2)  on  the  control  surface  or  end-plane 
and  LOS  that  enter  {a  > n/2)  on  the  control  surface  or  end  plane.  The 
conditions  that  determine  the  various  intersection  cases  are  listed  in 
Table  5.  When  the  intersection  case  has  been  determined,  the  coordin- 
ates of  the  entrance  (X^,  Z^)  and  exit  (X^,  Z^)  points  between  the  LOS 
and  the  plume  are  computed.  The  appropriate  assignments  are  indicated 
in  Table  5.  The  various  coordinate  definitions  are: 


s 
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Table  5.  LOS  - Plume  Intersection  Cases 


X = (L,  - z ) tan  a entrance  or  exit  on  end-plane 

e o 


X = -z  tan  at 
n o 


exit  on  nozzle-plane 


' 2 2 
o 

2 2 

R - R 
' m o 


entrance  or  exit  on  control  surface  (o'  ==  0,  tt) 


With  these  coordinates  defined,  the  total  length  of  the  path  through  the 
plume  is  computed  from 


(Xj  - + (Z,  - Z,) 


1 


The  total  length  of  the  path  within  the  plume  is  divided  into  N 
equal  size  intervals  and  an  array  of  path  position  defined  by 


s(i)  - (i  - 1)  i = 1.  N + 1 


Arrays  of  radial  and  axial  coordinates  corresponding  to  each  value  of  i 
are  computed  from 


r(i)  = + s(i)  sin  + d^ 


z(i)  = Z,  + s(i)  cos  or 


i = 1,  N + 1 


For  i - 1,  s(i)  = 0 and  r(i),  z(i)  are  the  cylindrical  coordinates  of  the 

entrance  point  of  the  LOS.  For  i - X t 1,  s(i)  = s and  r(i),  z(i)  are 

max 

the  coordinates  of  the  exit  point. 

Values  of  pTc  for  each  point  (i  - i,  X + 1)  along  the  path  are  com- 
puted by  linear  interpolations  based  on  the  input  plume  pTc  data  and  the 
coordinates  r(i),  z(i).  The  plume  data  f(n,  m)  [f  represents  p,  T,  or 
any  c]  are  defined  on  the  grid 

zz(n) n - 1,  NZ axial  static  is 

rr(n,  m)  ...  m = 1,  NR(n) radial  stations  at  each 

n =1,  NZ  axial  station 


The  first  step  of  the  interpolation  procedure  is  the  location  of  the  axial 
stations  z^  = zz(n),  - zz(n  + 1)  bounding  z(i),  i.  e.  , z^  < z(i)  < z^. 

The  procedure  now  diverges  depending  on  whether  r(l)  n or  r(i)  > 
where  is  the  maximum  grid  radius  at  z^,  r^^  = rrfn,  XR  (n)j.  For 
r(i)  <;  r^  (case  1),  the  radial  stations  at  z^  and  z^  that  bound  r(i)  are 
determined  and  intermediate  interpolations  performed  to  get  the  function 
f at  the  points  z^,  r(i)  and  z^,  r(i)  as  shown  in  Figure  ii-.  An  axial  inter- 
polation is  then  performed  using  f^,  f^,  z^,  z^,  and  z(i)  to  get  f at  z(i), 
r(i).  This  interpolation  procedure  can  be  interpreted  as  a construction 
of  the  f-surface  by  the  drawing  of  straight  lines  between  the  i profiles  at 
Zj  and  z^,  the  lines  being  equally  spaced  in  r out  to  r.  (sec  Figure  9). 
The  interpolation  procedure  for  r>  r^  (case  2)  is  based  on  the  comple- 
tion of  the  f- surface  construction  by  the  drawing  of  straight  lines  from 
the  point  r^  that  fan  out  to  the  f-profile  at  z^  (see  Figure  9).  The  inter- 
polation for  case  2 based  on  this  f-surface  construction  proceeds  by 
computing  the  radial  parameter  (see  Figure  10), 


r 

o 


r(i)  - i-j 
z(i)  - 
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and  performing  an  intermediate  interpolation  for  at  the  point  z^,  r^. 
An  axial  interpolation  is  then  made  for  the  final  value  of  f at  r(i),  z(i) 
by 


^ m 


whe  re 

f^  = f[n,  NR(n)]  = plume  boundary  value 

y = \/[z(i)  - + [r(i)  - rj^ 

^m  " ^^o  - ^1^^ 

If  the  computed  value  of  r^  is  greater  than  the  plume  radius  at  z^,  i.  e.  , 
= rr[n  + 1,  NR(n  + 1)],  then  the  point  z{i),  r(i)  lies  in  the  shaded  area 
shown  in  Figure  10.  Although  the  point  is  inside  the  plume  control  vol- 
ume, the  value  of  f is  set  to  the  plume  boundary  value,  f^. 
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E.  Subroutine  DELTRAN 


Subroutine  DELTRAN  is  the  primary  band  model  calculation  rou- 
tine. The  input  consists  of  arrays  of  pTc  data  and  band  model  para- 
meters describing  an  optical  path.  Specifically,  the  following  arrays 
are  entered 


s(i) 

P(i) 

c(i) 

k(i) 

D(i)  = l/6(i) 
Yt  (i) 


path  position  (cm),  s(i)  = 0 
pressure  (atm) 
concentration  (mole  fraction) 
absorption  parameter  (cm  Vatm) 
line  density  (1/cm  ^) 

Lorentz  line  width  parameter  (cm  ^) 
Doppler  line  width  parameter  (cm 


>-  i = 1,  N + 1 
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Mode  data  are  also  entered  to  indicate  line  shape  and  inhomogeneity 
approximation  to  be  used.  The  return  of  the  subroutine  are  the  trans- 
mittance T(i)  and  transmittance  derivative  dT(i)/ds  for  a single  species. 
The  method  of  computation  has  been  given  in  Section  II.  Path  integrals 
appearing  in  the  band  model  formulation  are  computed  by  trapezoidal 
quadrature  on  the  coordinate  axis  s{i). 

For  a given  species,  three  calls  to  DELTRAN  are  made  by  the 
main  program  ATLES.  The  first  call  computes  '(i)  and  dT(i)/ds  along 
the  absorption  path.  For  these  calculations  the  initial  values  of  path 
integrals  are  set  to  zero  at  i = 1.  When  the  calculation  loop  ends  at 
i = N + I,  the  cumulative  sums  representing  the  path  integrals  over 
the  entire  absorption  path  are  saved.  The  second  call  to  DELTRAN 
computes  T(i)  and  dT(i)/ds  over  the  emission  portion  of  the  total  optical 
path  as  though  the  absorption  path  were  not  present.  This  calculation 
is  carried  out  by  again  initiating  integral  sums  to  zero  at  i = 1.  The 
cumulative  sums  for  this  calculation  are  not  saved.  The  third  call  to 
DELTRAN  computes  T(i)  and  dx(i)/ds  over  the  emission  path,  but  this 
time  with  account  of  absorption  by  the  coupled  absorption  path.  This 
calculation  is  performed  by  initializing  the  path  integrals  at  i = 1 to  the 
cumulative  sums  that  have  been  saved  from  the  first  call  to  DELTRAN. 


IV.  PREPARATION  OF  INPUT  DATA 


A computational  run  of  ATLES  requires  a set  of  program  control 
cards  (Figure  1 1 ) to  specify  the  mode  of  computation  and  to  supply  input 
data.  A control  card  may  also  signal  that  auxiliary  data  are  to  be  read 
in  by  the  program.  These  auxiliary  data  may  be  provided  in  card-deck 
form  or  as  attached  tape  or  disc  files.  When  the  set  of  control  cards 
and  auxiliary  data  are  read  in,  a computational  run  is  made  and  an  output 
listing  is  generated.  The  program  then  continues  to  read  a new  set  of 
program  control  cards  and  execute  a new  run.  The  computations  of  the 
second  and  subsequent  runs  proceed  with  the  conditions  of  the  first  run 
as  modified  by  changes  indicated  by  the  subsequent  control  cards. 

Each  type  of  control  card  contains  an  alphanumeric  name  in  the 
first  ten  card  columns,  and  must  be  left-justified.  If  data  are  specified 
on  the  control  card,  they  must  be  entered  in  accord  with  the  format 
specification  indicated  in  the  detailed  descriptions  given  below.  All 
fields  of  the  program  control  cards  are  10  columns  wide.  In  general, 
integer  and  alphanumeric  (Hollerith)  data  must  be  right-justified  in  their 
fields.  Noninteger  numerical  data  may  be  entered  in  either  F-format 
with  decimal  point  or  right-justified  E-format  with  decimal  point  and 
exponential  symbol  E.  The  types  of  control  cards  and  the  data  contained 
on  them  are  illustrated  in  Figure  11.  Detailed  descriptions  of  each  card 
type  are  given  in  the  following  sections, 

A.  Title  Card;  The  card  name  is  TITLE.  Columns  11-60  oi  this  card 
may  be  used  for  any  identification  title  desired.  The  title  will  be 
displayed  in  the  final  listing  of  the  computation  run  and  on  any  plots 
generated  by  the  run. 

B.  Calculation  Mode  Card;  The  card  name  is  CAl.CMODE.  The  varia- 
ble CMODE  (format  AlO)  must  have  one  of  the  three  alphanumeric 
values  ATM,  SRC,  or  CPL.  The  value  ATM  signals  that  calculation 
for  an  absorption  path  only  will  be  made.  This  mode  would  be  used, 
for  example,  in  computing  the  transmittance  spectrum  over  an 
atmospheric  slant  path.  Similarly,  the  value  SRC  indicates  that 
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calculations  for  an  emission  path  only  will  be  made.  If  the  value 
is  CPL,  full  calculations  for  a coupled  emission/absorption  optical 
path  are  performed.  If  the  emission  source  is  specified  as  a plume 
(see  SOURCE  card)  and  the  absorption  path  as  an  atmospheric  slant 
path  (see  ATMOSPHERE  card),  for  example,  calculations  would  be 
made  for; 

1.  Plume  radiance  at  plume  boundary 

2.  Plume  radiance  attenuated  by  atmosphere 

3.  Transmittance  (and  radiance)  of  atmospheric  slant  path 

4.  Effective  atmospheric  slant  path  transmittance 

The  variable  LSHAPE  (format  A 10)  on  the  CALCMODE  card 
defines  the  spectral  line  shape  that  will  be  used  in  the  statistical 
band  model  calculations.  This  variable  must  have  one  of  the  three 
alphanumeric  values  LORENTZ,  DOPPLER,  or  VOIGT.  For  cal- 
culations involving  atmospheric  slant  paths,  the  use  of  the  Lorentz 
line  profile  is  generally  adequate  if  one  end  of  the  slant  path  is 
below  20km. 

Path  Intervals  Card;  The  card  name  is  INTERVALS.  The  variables 
NAPATH  and  NSPA  IH  (each  format  110)  are  the  number  of  equal 
geometric  length  intervals  into  which  the  absorption  and  emission 
path,  respectively,  are  divided  for  numerical  integrations  over  the 
paths.  The  minimum  allowed  value  of  either  variable  is  one  and 
the  maximum  allowed  value  is  100.  The  value  of  these  two  variables 
must  be  large  enough  to  give  reasonable  numerical  accuracy,  but, 
in  order  to  minimize  computer  cost,  should  not  be  too  much  larger 
than  necessary.  For  a large  run  in  which  calculations  will  be  made 
for  many  lines  of  sight  through  a plume  for  many  spectral  positions, 
an  empirical  determination  of  NAPATH  and  NSPATH  can  be  made 
for  one  line  of  sight  at  one  spectral  position  (see  discussion  of 
SPECTRUM  and  GRID  program  control  cards).  The  spectral  posi- 
tion selected  should  be  one  for  strong  absorption  and  the  line  of 
sight  should  be  one  which  passes  through  the  plume  where  the  opti- 
cal density  is  large.  Then,  by  making  a series  of  runs  with 

Except,  an  equal  mass  division  is  made  for  atmospheric  slant  paths 


I 


Increasing  NAPATH  and  NSPATH,  optimum  values  can  be  deter- 
mined by  noting  how  the  computed  results  converge. 


For  problems  involving  a plume  source,  it  may  not  be 
desirable  to  use  a fixed  value  for  NSPATH.  A fixed  value  implies 
that  the  portion  of  the  line  of  sight  through  the  plume  is  divided 
into  the  same  number  of  integration  intervals  w'hether  the  line  of 
sight  passes  through  a full  diameter  of  the  plume  or  just  grazes 
the  boundary.  In  the  former  case,  a large  value  of  NSPATH  may 
be  required  whereas  in  the  latter,  a value  near  one  would  probably 
be  sufficient.  An  effective  variable  value  of  NSPATH  can  be  used 
for  such  plume  calculations  by  setting  NSPATH  to  zero  (or  blank) 
and  specifying  the  approximate  source  path  increment  directly. 

This  increment  (unit  = m)  is  entered  as  DELS  (format  ElO)  on  the 
interval  card.  Integrations  along  the  source  portion  of  the  line  of 
sight  will  be  made  with  an  integration  interval  approximately  equal 
to  DELS  regardless  of  the  length  of  the  optical  path  through  the 
plume. 

For  calculations  involving  absorption  or  emission  cells  (see 
SOURCE  and  ATMOSPHERE  cards)  NAPATH  and  NSPATH  may  be- 
set (appropriately)  to  zero.  In  this  case,  equal  divisions  along  the 
optical  path  segments  are  not  made,  but  rather  the  coordinates 
used  to  enter  the  pTc  data  along  the  path  segments  are  used  to  de- 
fine the  integration  interval  subdivision. 

D.  Inhomogeneity  Approximations  Card:  The  card  name  is  INHOMMODE. 

The  variables  MODEATM,  MODESRC,  and  MODECPL  specify  the 
approximation  that  will  be  used  to  handle  nonuniformities  along  the 
absorption,  emission,  and  coupled  absorption/ emis s ion  paths, 
respectively.  The  format  of  each  is  AlO,  and  the  allowed  values 
are  UNIFORM,  CG  and  DR.  The  value  CG  indicates  that  the 
Curtis-Godson  approximation  will  be  used  while  DR  indicates  that 
the  derivative  approximation  will  be  used.  MODEATM  = CG  is 
adequate  for  atmospheric  slant  paths  and  may  be  adequate  for 
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computing  source  radiance  If  thermal  gradients  are  not  too  large. 
For  large  thermal  gradients  within  the  source,  MODESRC  should 
be  set  to  DR.  Computation  of  hot  gas  emission  attenuated  by  an 
intervening  atmospheric  path  should  almost  always  be  made  with 
MODECPL  = DR. 

The  value  UNIFORM  indicates  that  calculations  will  be  made 
as  though  the  path  were  uniform,  and  is  an  appropriate  value  only 
for  MODEATM  and  MODESRC.  If  MODEATM  and/or  MODESRC 
have  this  value,  an  override  of  the  values  NAPATH  and/or  NSPATH 
on  the  INTERVALS  card  is  made  and  calculations  arc  made  as 
though  they  had  been  specified  as  one.  The  entered  value  of 
NSPATH  still  gives  the  subdivision  of  the  source  path  segment  used 
for  the  coupled-path  calculation. 

E.  Spectrum  Card:  The  card  name  is  SPECTRUM.  The  variable 

UNIT  (format  AlO)  determines  the  spectral  units  required  for  input 

spectral  data  and  also  the  spectral  unit  of  computed  results.  If 

UNIT  has  the  value  MICRON,  input  spectral  data  must  be  in  um, 

2 

and  computed  radiance  - suits  will  be  expressed  in  W / cm  -sr-pm 

as  a function  of  wavelength  in  pm.  If  UNIT  has  any  other  value. 

spectral  data  must  be  entered  as  wavenumber  with  unit  cm 

2 - 1 

Computed  radiance  spectra  will  be  listed  with  unit  W/cm  -sr-cin 
as  a function  of  wavenumber  with  unit  cm 


The  variables  WMIN  and  WMAX  are  the  lower  and  upper 
extent,  respectiv'ely,  of  the  spectral  region  to  be  considered. 

The  format  specification  of  each  is  ElO  and  the  unit  must  corres- 
pond to  UNIT.  The  region  between  WMIN  and  WMAX  is  divided 
into  NSI’EC  (format  110)  equal-size  intervals  and  calculations  per- 
formed at  the  resulting  NSPEC  + 1 spectral  positions.  The  mini- 
mum and  maximum  allowed  values  for  NSPEC  arc  0 and  200, 
respectively.  If  NSPEC  is  0,  calculations  at  only  the  one  spectral 
position  WMIN  are  made. 
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Bandpass  Card;  The  c ard  name  is  J^ANDPASS.  The  variables  WL 
and  WU  (each  format  EiO)  are  the  lower  and  upper  extents,  res- 
pectively, of  a spectral  bandpass  over  which  the  computed  radiance 
and  transmittance  spectra  will  be  a\-e raged.  In  order  to  be  mean- 
ingful, these  limits  should  be  the  same,  or  inside,  the  spectral 
limits  set  by  the  SPECTRUM  card.  The  BANDPASS  card  can  be 
used  to  compute  average  spectral  values  within  a specified  sensor 
bandpass  if  the  inband  responsivity  of  the  sensor  is  assumed  to  be 
constant.  Up  to  10  different  BANDFUASS  cards  may  he  vtsed  in  one 
computational  run.  The  unit  of  VVL  and  \VU  must  be  consistent  with 
the  value  of  UNIT  on  the  SPECTRUM  card. 

The  BANDPASS  card  may  also  be  used  to  eliminate  all  pre- 
viously defined  bandpasses  by  giving  WL  the  alphanumeric  value 
RESET. 

G.  Geometry  Card;  The  card  name  is  GEOMETRY  and  is  used  to 
define  the  geometry  of  observation  along  an  atmospheric  slant 
path.  This  geometry  may  be  defined  in  one  of  two  ways. 

If  PCODE  (format  A 10)  has  the  value  EXOATM,  the  geometry 
consists  of  a sensor  above  the  atmosphere  looking  down  on  a target 
plume  within  the  atmosphere  (sec  Figure  3).  In  this  case,  the 
variable  ALT  (format  EIO)  is  the  altitude  of  the  targei.  plum.e  abo\e 
ground  level  (unit  = km),  and  ZENITH  (format  EIO)  is  the  zenith 
angle  (unit  = degrees)  from  the  plume  to  the  sensor  (ZENITH  = 0 
implies  that  the  sensor  is  directly  above  and  looking  dovt'n  on  tlie 
target). 

If  F’CODE  has  the  value  ENDOAFM,  the  geometry  is  defined 
with  reference  to  a sensor  within  the  atmosphere  (see  Figure  4). 
The  omission  source  may  be  cither  within  or  abo\e  the  atmosphere. 
In  this  case,  ALT  is  the  sensor  altitude  and  ZENITH  is  the  zenith 
angle  from  the  sensor  to  the  source  (ZENITH  -■  0 implies  that  the 
sensor  is  directly  below  and  looking  up  at  the  source).  The 
variable  RANGE  (format  EIO)  is  tin-  range  in  km  from  the  sensor 
to  the  target  and  needs  to  he  defined  only  if  lYiODE  - F.'NDOATM. 
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The  variable  ZMAX  (format  ElO)  defines  the  "top  of  the 
atmosphere"  and  may  be  any  value  less  than  or  equal  to  100km. 

The  unit  is  km. 

If  the  emission  source  is  a plume  (see  SOURCE  card),  the 
variable  ASPECT  (format  ElO)  is  used  and  is  the  viewing  angle 
between  the  plume  axis  and  the  line  of  sight.  The  unit  is  degrees, 
and  the  angle  is  measured  from  nose-on  to  the  missile  or  aircraft 
generating  the  plume. 

H.  Print  Card:  The  card  name  is  PRINT  and  contains  the  two  variables 

WATM  and  WSRC  (each  format  AlO).  If  WATM  has  the  value 
PRINT,  a detailed  printout  of  all  major  computational  steps  through- 
out the  program  for  the  absorption  path  portion  of  the  total  optical 
path  results.  Similarly,  if  WSRC  = PRINT,  printouts  for  the  source 
portion  occur.  The  use  of  this  card  should  be  restricted  to  debug- 
ging or  computational  check  runs  only.  If  any  more  than  one 
spectral  position  for  one  line  of  sight  is  involved  in  the  job  run, 

the  printout  will  be  enormous. 

I.  Plot  Cards  ; The  two  cards  named  PLOTRAD  and  PLOTTRAN 
specify  if  any  film  plots  are  to  be  made  of  the  computed  radiance 
and  transmittance  spectra,  respectively.  The  variable  WSCALE 
(format  ElO)  is  the  spectral  plot  scaling  parameter  and  has  unit 
cm  ^/(inch  of  plot)  or  g,m/(inch  of  plot)  depending  on  the  value  of 
UNIT  on  the  SPECTRUM  card.  The  value  of  WSCALE  should  be 
chosen  so  that  the  entire  spectral  interval  specified  by  the 
SPECTRUM  card  is  contained  within  a 12-inch  plot.  Note  that 
WSCALE  is  included  on  both  the  radiance  and  transmittance  plot 
cards.  The  last  value  of  WSCALE  read  is  the  value  used. 


This  description  is  of  concern  only  to  users  at  The  Aerospace  Corpora- 
tion since  the  plot  functions  used  in  the  coding  are  in-house  routines. 
Users  outside  The  Aerospace  Corporation  must  write  their  own  plot 
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If  a PLOTRAD  card  is  included  in  the  input,  but  CSRC, 

CSNSR  and  CATM  are  not  specified  (normal  usage),  the  following 
radiance  plots  will  be  generated. 

1.  Source  radiance  with  an  optimum  y-scale  factor 

2.  Sensor  radiance  (attenuated  source  radiance)  wi.th  an  optimum 
y-scale  factor 

3.  Sensor  radiance  with  the  same  y-scale  factor  as  plot  one 
(but  only  if  the  factor  is  different  from  that  for  plot  two) 

4.  Atmospheric  radiance  with  an  optimurri  y-scale  factor 

5.  Atmospheric  radiance  with  the  same  y-scale  factor  as  plot 
one  (but  only  if  the  factor  is  different  from  that  for  plot  four) 

If  CSRC,  CSNSR  or  CATM  (each  format  ElO)  are  specified,  addi- 
tional plots  for  the  source,  sensor,  or  atmospheric  radiance, 
respectively,  are  generated  using  these  values  as  the  y-scaling 
factors.  The  unit  of  these  variables  is  (spectral  radiance  consis- 
tent with  UNIT)/(inch  of  plot).  If  CSRC,  CSNSR  or  CATM  are  the 
same  as  the  scale  factor  for  any  of  the  five  basic  radiance  plots, 
no  additional  plots  are  made.  The  maximum  allowed  plot  height  is 
10  inches. 

The  PLOTTRAN  card  operates  in  a similar  manner.  The 
basic  transmittance  plots  are 

1.  Atmospheric  transmittance  with  a y-scale  factor  of  0.  10 

2.  Effective  atmospheric  transmittance  with  a y-scale  factor  of 

0.  10 

If  CT  and  CTE  are  not  zero,  additional  transmittance  plots  arc 
made  for  the  atmospheric  and  effective  atmospheric  t ransmittances , 
respectively,  using  these  values  as  the  y-scale  factor.  The  unit  of 
CT  or  CTE  is  (full  scale  percent  t ransmi ttanc e / 1 000 )/( inch  of  plot) 
for  a 10  inch  high  plot. 
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J.  Species  Card;  The  card  name  is  SPECIES  and  is  used  to  include 
a particular  molecular  species  in  the  calculations  and  to  call  for 
the  input  of  band  model  parameters  for  that  species.  SNAME 
(format  AlO)  is  an  arbitrary  identification  name  for  the  species 
(e.g.,  H20,  C02,  N20)  and  SNUMBER  (format  110)  is  an  identifi- 
cation number  (allowed  values  are  1,  2 or  3).  If  SNAME  has  the 
value  DELETE,  the  species  identified  by  SNUMBER  is  deleted 
from  current  and  subsequent  calculations.  If  NFILE  is  zero  or 
blank,  band  model  parameters  for  the  species  must  follow  the 
SPECIES  card  immediately  in  the  input  card  file.  The  required 
deck  structure  is  given  in  Figure  12.  If  NFILE  is  not  zero,  then 
the  value  specifies  the  section  on  the  Band  Model  Parameter  File 
containing  the  desired  parameters.  This  file  must  have  been 
attached  to  the  job  run  on  input  device  labeled  2 ' . The  band  model 
parameter  section  identification  numbers  are  given  in  Table  2.  If 
the  variable  PRINT  (AlO)  has  the  value  PRINT,  a listing  of  the 
band  model  parameters  will  be  made.  Up  to  three  SPECIES  cards 
(with  different  SNUMBER)  may  be  used  in  one  computational  run. 

K.  Atmosphere  Card;  The  card  name  is  ATMOSPHERE.  The  varia- 
ble ACODE  (format  AlO)  must  be  one  of  the  three  values  ATM,  ■* 

LOS  or  NOREAD.  The  value  ATM  indicates  that  the  absorption  ^ 

path  is  an  atmospheric  slant  path  and  that  data  specifying  the  alti- 
tude profile  of  a model  atmosphere  should  now  be  made  available. 

If  NFILE  (format  110)  is  blank  or  zero,  the  model  atmosphere  data 
is  assumed  to  be  in  card  deck  form  and  to  immediately  follow  the 
ATMOSPHERE  control  card.  The  required  deck  structure  is  given 
in  Figure  13.  If  NFILE  is  not  zero,  the  value  indicates  the  section 
of  the  Model  Atmosphere  File  that  contains  the  desired  model.  This 
file  must  have  been  previously  attached  to  the  Job  on  input  de\  ice  3 . 


For  the  CDC  7600  computer  with  SCODE  3 operating  system,  the*  Band 
Model  Parameter  File  must  have  been  attached  to  the  job  with  file  nanu- 
TAPE2  and  the  Model  Atmosphere  File  with  file  name  TAPES. 
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NAME 

Identification  name  (format  A10)=i= 

Av 

Spectral  resolution  of  data  (cm'^)^ 

ad) 

Ratio  of  resonant  self-broadening  parameter  to 

Y 

a(2)t 

Ratio  of  H^O  foreign  gas  broadening  parameter 

to  Y 

a(3)t 

Ratio  of  CO2  foreign  gas  broadening  parameter 

to  Y 

a(4) 

Ratio  of  N^  foreign  gas  broadening  parameter  to  y 

a(5) 

Ratio  of  ©2  foreign  gas  broadening  parameter  to  y 

a(6) 

Mass  of  absorbing  species  (amu) 

NW 

Number  of  elements  in  wavenumber  array  (format  15) 

NT 

Number  of  elements  in  temperature  array  (forn 

.at  15) 

v(i) 

Wavenumber  array  (cm ’^)i=  1,2,3,  ...,NW 

T(j) 

Temperature  array  (°K)  j = 1,2,3,...,  NT 

I^(i,  j) 

Mean  absorption  coefficient  array  (cm'Vatm) 

D(i.  j)§ 

Effective  line  density  paranieter  array  (1/cm  ^ 

) 

Y(‘) 

Broadening  parameter  for  non-resonant  self- 
broadening at  STP  (cm"Vatm) 

Unless  otherwise  indicated,  as  here,  all  fields  are  format  ElO  and  data 
contain  appropriate  decimal  points  and  exponential  symbols  E. 

^The  interval  size  between  elements  of  the  v(i)  array  may  be  Av  also, 
but  need  not  be. 

^If  the  parameter  set  is  for  H2O,  a(2)  =0.  If  for  CO2.  - 0. 

SD(i,  j)  = 1/6  (i.  j).  This  parameter  is  the  effective  line  density  parameter 
consistent  for  use  in  the  statistical  band  model  for  equal  strength  lines. 


Figure  12.  Band  Model  Parameter  Set  Deck  Structure 


-6  5- 


I 


IDNAME 

model  atmosphere  identification  name 
(NAME  is  optional,  but  card  must  be  included) 

h(i) 

altitudes  (m).  h(i)  must  be  the  33  values 
0, 1,  2 24,  25,  30,  35, 40,  45,  50,  70,  100  km. 

T(i) 

temperature  (K) 

lUNIT 

units  flag 

If  lUNIT  = 1 

p(i)  - pressure  (atm) 

Cj^  ~ water  vapor  concentration 

(mole  fraction) 

If  lUNIT  1 

p(i)  - pressure  (mb)  ^ 

Cr  q(>-)  - water  vapor  concentration  (g/m  ) 

Figure  13.  Deck  Structure  for  Model  Atmosphere 
Input  Data. 
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The  section  identification  numbers  are  given  in  Table  1 . If  the 
variable  PRINT  (format  AlO)  has  the  value  PRINT,  the  model 
atmosphere  data  will  be  listed.  The  variables  Cl,  C2  and  C3  are 
used  to  specify  concentration  conditions  for  species  identified  as 
1,  2 and  3,  respectively,  by  the  SPECIES  cards.  If  Ci  is  zero  or 
blank  (format  ElO),  species  i is  not  considered  ,-^s  an  atmospheric 
constituent.  If  Ci  has  a numerical  value  (format  ElO)  this  value 
will  be  used  as  the  constant  mixing  ratio  of  species  i in  the  atmos- 
phere, If  Ci  has  the  value  H20  (format  AlO),  the  water  vapor  con- 
centration profile  of  the  model  atmosphere  will  be  used  to  determine 
concentrations  for  species  i. 

A user  may  specify  the  pressure,  tem.perature  and  species 
concentration  (pTc)  along  an  absorption  line  of  sight  by  setting 
ACODE  to  the  value  LOS.  In  this  case,  NTAPE  should  be  set  to 
five  and  a card  deck  of  the  pTc  data  inserted  immediately  after  the 
ATMOSPHERE  card.  The  required  input  card  deck  structure  is 
given  in  Figure  14.  If  PRINT  has  the  value  print,  this  pTc  data 
will  be  listed.  The  pTc  deck  allows  for  the  specification  of  up  to 
five  species.  The  variables  Ci  (format  110)  in  this  case  identify 
the  "species  column"  containing  the  concentration  data  to  be 
assigned  to  species  i.  "Species  cohunn"  is  defined  by:  species 

column  1 = card  columns  31-40,  species  column  2 - card  columns 
41  -50,  . . . , species  column  5 = card  columns  71  -80.  For  example, 
if  species  2 has  been  identified  by  the  SPECIES  card  to  be  CO,  and 
the  CO  concentration  has  been  entered  on  the  pTc  data  deck  in  card 
columns  61-70,  then  C2  = 4. 


^7' 


'4^: 

.y 


The  value  ACODE  = NOREAD,  indicates  that  neither  an 
atmospheric  model  nor  a pTt  data  deck  is  to  be  read  in  (presumably 
one  of  these  having  been  read  in  on  a previous  run),  bvit  that  new 
data  assignments  ar<-  tc>  be  macle  by  a new  set  of  Ci  variables. 

Source  Card:  The  card  name  is  SOURCE  and  operates  in  much  the 

same  way  as  the  ATMOSPHERE  card  does.  SCODE  (format  AlO) 
may  be  one  of  the  three  values  PLM,  LOS  or  NOREAD.  The  first 
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IDNAME 

N 

NAMEl  thru 
NAME5 

s(i) 

PU) 

T(i) 


cell  identification  name  (optional) 

number  of  path  positions  for  which  data  is  to  be 
entered.  2 i N £ 101  and  N is  integer 

species  identification  names  (optional) 


path  positions.  s(l)  - 0 and  s(i)  increases  with  i 
pressure  (atm) 
temperature  (K) 

species  concentrations  (mole  fraction). 

1 < m £ 5 


Figure  14.  Deck  Structure  for  Emission  or  Absorption  Cell  Input  Data 
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value  IndLcates  that  the  radiation  source  to  be  considered  is  a 

plume  and  the  second  that  the  source  is  a single  line  of  sight 

through  a hot  gas  region.  The  nominal  input  mode  for  the  plume 

or  line  of  sight  pTc  data  is  a card  deck  immediately  following  the 

SOURCE  card.  The  structures  for  the  line  of  slight  and  plume  ! 

pTc  data  decks  are  given  in  Figures  14  and  15,  respectively.  For 

card  input,  NTAPE  (format  110)  must  be  5.  If  the  pTc  data  J 

are  attached  to  the  job  in  tape  or  permanent  file  form,  they  should  | 

be  attached  on  input  device  i and  NTAPE  set  to  i.  A value  i = 7,  8 

or  9 is  allowed.  If  PRINT  (format  A 10)  has  the  value  PRINT,  the 

pTc  input  data  will  be  listed.  The  parameters  Cl,  C2  and  C3  (each 

format  110)  are  the  species  column  in  which  the  concentration  data  .■ 

for  species  1,  2 and  3,  respectively,  are  entered  in  the  pT c data  i 

deck.  The  value  SCODE  = NOREAD  is  used  when  only  the  species  i 

column  assignment  is  to  be  changed. 

M.  Plume  Integration  Card;  The  CARD  name  is  GRID  and  is  used  to 
specify  how  the  Integration  over  the  projected  area  of  a plume  is 
to  be  performed.  Two  forms  of  the  card  are  allowed.  If  calcula- 
tions for  a single  line  of  sight  through  a plume  are  desired,  the  ' 

NTAPE  field  (format  110)  should  be  left  blank.  The  line  of  sight 
is  then  defined  by  the  axial  distance  Z (format  ElO  and.  unit  = m)  / 

from  the  nozzle  and  the  transverse  distance  D (format  ElO  and  } 

unit  = m)  from  the  plume  axis  of  the  point  of  intersection  between 
the  line  of  sight  and  the  projected  area  of  the  plume.  If  NTAPE  = i * 

(format  110)  is  not  blank,  the  Z and  D coordinates  defining  many 
lines  of  sight  through  the  plume  will  be  read  in  from  input  device  i. 

If  NTAPE  - 5,  the  grid  data  must  be  in  card  deck  form  and  immed- 
iately follow  the  GRID  card.  The  required  deck  structure  is  given 
in  Figure  16.  If  the  variable  PRINT  (forrr^at  A 10)  has  the  \alue 

PRINT,  the  grid  data  will  be  listed.  If  the  \ariable  WMAP  (format  i 

j 

AlO)  has  the  value  PRINT,  the  computational  rc^sults  of  the  bandpass  , 

averaged  radiance  and  transmittance  for  each  line  of  sight  through 
the  plume  will  be  listed.  (The  normal  output  is  the  area  averaged 
result  for  all  the  lines  of  sight.  ) 

J 

i 


-69- 


I 


i 


\ 


I 


¥ 
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IDNAME 


plume  identification  name  (optional) 


L 


control  volume  length  (m) 


R nozzle  radius  (m) 

o 

R maximum  control  volume  radius  (m) 

m 


NZ 


number  of  axial  stations  for  which  data  is  to  be 
entered,  2 £ NZ  5 31  and  NZ  is  integer 


NAMEl  thru  NAMES  species  identification  names  (optional) 


Z(i) 


NR(i) 


r(i,j) 


P(i.  j) 


I(Lj) 


axial  station  positions  (m). 

Note:  Z(l)  = 0 (exit  plane),  Z(NZ)  = L,  and 
Z(i)  increases  with  i 

number  of  radial  stations  at  axial  station  i. 

2 < NR(i)  <11  and  NR(i)  is  integer 

radial  station  positions  (m).  Note:  r(i,  1)  = 0 (axis). 
r[i.  NR(i)J  - - R^)  + R^,  r(i.j)  increases 

with  j 

pressure  (atm).  If  p(l,  1)  = 0,  the  ambient  atmo- 
spheric pressure  at  the  plume  altitude  is  used  for 

all  p(i,j) 

temperature  (K) 

species  concentrations  (mole  fraction).  1 < ni  < 5 


Note:  The  boundary  values  of  p(i,j),  T(i,j)  and  c (i,j)[i.e.,  for  j NR(i)] 

should  be  set  to  the  ambient  atmospheric  va\ues. 


Figure  15.  Deck  Structure  for  Plume  Input  Data 
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IDNAME 

grid  identification  name  (optional) 

NZ 

number  of  axial  grid  coordinates  where  LOS 
calculations  will  be  made.  1 S NZ  s 21  and  NZ 

is  integer 

/(!)'■' 

axial  grid  coordinate  (m).  Z(i)  increases  with  i 

ND(i) 

number  of  transverse  grid  coordinates  at  axial 
coordinate  i.  1 i ND(i)  5 8 

n(i,j) 

transverse  grid  coordinates  (m).  D(i.  j)  in- 
creases with  j.  If  ND(i)  £ 6,  delete  second  card 

This  coordinate  is  measured  in  the  plume  coordinate  systeni,  the 
projected  area  coordinate  system.  For  aspect  angles  near  0 or  I80“. 
negative  values  of  / may  be  required  in  order  to  properly  specify  the 
desired  line  of  sight. 


Figure  16.  Deck  Structure  for  Integration 
Grid  Input  Data 
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The  following  integration  cases  can  be  handled  by  a suitable 
combination  of  the  GRID  control  card  and  auxiliary  grid  data. 

I 

Case  1.  Single  Line  of  Sight 

In  this  case,  only  the  GRID  card  (Type  1)  is  required.  The 
coordinates  Z(l)  and  D(l,  I)  are  obtained  directly  from  the  GRID 
card.  NZ  and  ND(I)  are  automatically  set  to  1.  If  ASPECT  = 0 or 
180  (see  GEOMETRY  card),  Z(l)  need  not  be  specified. 

Case  2.  Single  Axial-Station 

The  auxiliary  grid  data  deck  consists  of  one  or  two  card  with  XZ 
specified  as  1 and  containing  the  D(l,j)  values  where  computations 
are  to  be  made.  Note:  D(l,  1)  must  be  zero. 

Case  3.  Whole  Plume  Average  (ASPECT  4 ^ or  180) 

The  full  auxiliary  data  deck  is  required.  Note:  D(i,  1)  must 
be  zero  for  all  i. 

Case  4,  Whole  Plume  Average  (ASPECT  = 0 or  180) 

This  case  is  the  same  as  Case  2 except  that  Z(l)  need  not  be 
specified,  and  D(l,  1)  need  not  be  zero.  For  example,  if  ASPECT 
= 0,  D(l,  1)  should  be  set  equal  to  the  obscuration  disc  radius  R^. 

For  an  integration  over  the  whole  projected  area  of  a plume, 
prior  calculations  are  recommended.  Although  the  plume  is  defined 
throughout  a paraboloid  control  volume,  the  area  of  integration 
specified  by  the  grid  data  need  not  (and  for  economy's  sake,  should 
not)  be  the  entire  projected  area  c''  this  volume.  The  integration 
area  need  be  only  the  region  where  the  majority  of  the  radiation 
is  emitted.  A series  of  runs  using  the  Type  1 GRID  program  con- 
trol card  with  D(l,  1)  - 0 can  be  used  to  estimate  the  axial  extent  of 
the  integration  area.  These  computations  can  usually  be  done  for 
one  wavenumber.  The  transverse  extent  of  the  significant  part  of 
the  plume  at  each  axial  station  can  be  estimated  from  the  input 
plume  data  profiles.  Or,  radiance-  profiles  at  a selected  axial 
stations  can  be  generated  by  a series  of  runs  with  the  Type  1 GRID 
card  with  constant  Z and  increasing  D. 
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When  the  desired  extent  of  the  plume  integration  area  has 
been  determined,  a full  area  integration  run  can  be  done.  The 
integration  over  the  projected  area  of  che  plume  uses  a simple 
trapezoidal  numerical  quadrature  routine.  The  grid  should  be 
structured  with  a spatial  resolution  fine  enough  to  justify  this 
approach. 

The  specification  of  how  the  spatial  integration  over  the 
projected  area  of  a plume  should  be  carried  out  is  the  most  tedious 
problem  in  data  preparation,  particularly  for  viewing  aspect  angles 
which  are  not  either  end-on  or  broadside.  The  basic  problem  is  to 
choose  a grid  of  points  on  the  projected  area  of  the  plume  that 
adequately  covers  the  areas  of  dominant  radiation,  represents  a 
fine  enough  mesh  so  that  the  numerical  integration  is  accurate, 
and  is  limited  enough  so  that  the  cost  of  computation  is  not  exces- 
sive. Attempts  to  automate  the  selection  of  an  integration  grid 
have  been  made  but  the  routines  are  not  yet  considered  useful  enough 
to  be  incorporated  into  the  code.  Thus,  the  specification  of  the 
integration  grid  is  left  to  the  user's  discretion.  An  example  grid 
selection  is  given  in  Section  V. 

N.  Execution  Card:  The  card  name  is  RUN.  When  a RUN  card  is 

encountered,  computations  are  begun  using  the  data  entered  up  to 
that  point,  and  an  output  listing  of  the  results  is  made.  When  the 
computation  and  listing  are  completed,  the  program  continues  to 
read  program  control  cards  until  a new  RUN  card  is  encountered. 

A new  computation  is  then  begun  using  all  of  the  conditions  and 
data  of  the  previous  run  except  those  which  have  been  changed,  added 
or  deleted  by  the  intervening  program  control  cards  and  auxiliary 
data  decks.  The  process  is  repeated  until  an  end-of-file  card  is 
encountered.  With  this  feature,  several  different  computational 
runs  can  be  made  with  one  job  submission.  Other  than  the  necessity 
that  all  required  data  be  entered  before  a RUN  card  is  encountered 
and  that  auxiliary  data  in  card  form  immediately  follow  the  con- 
trol card  that  calls  for  them,  the  program  control  cards  niay  be 
f in  any  ordc r . 
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Great  care  should  be  taken  In  the  preparation  of  input  data 
since  there  are  very  few  internal  checks  of  data  consistency.  A 
general  feature  in  data  preparation  is  that,  if  particular  data  on  a 
control  card  i re  not  required,  they  need  not  be  specified.  If  none 
of  the  data  on  a control  card  are  required,  that  card  nee  a not  be 
included.  > 


V.  EXAMPLE  CALCULATIONS 


A.  ERIM  Hot-through  Cold  H2O  Spectra 

An  experimental  study  designed  specifically  to  measure  the 

radiance  spectra  of  hot  gases  as  viewed  through  cool  absorption  paths 

(12) 

has  recently  been  completed  at  ERIM.  These  measurements  pro- 

vide quantitative  data  against  which  band  model  formulations  can  be 
tested.  This  example  calculation  illustrates  the  use  of  ATLES  in  pre- 
dicting radiance  spectra  for  experimental  set-ups  such  as  used  by 
ERIM, 

The  ERIM  runs  considered  are  identified  as  9R  and  IIMR,  The 
hot  emission  source  is  a hot-gas  cell  containing  H^C  and  N^-  The 
optical  path  through  the  source  is  uniform  with  L = 60cm,  T = 1200°K, 
p = 0.  10  atm,  and  c(H^O)  = 0.  50.  The  radiation  from  the  hot  H^O  in 
the  2.  7-pm  band  is  viewed  through  a long  cool  absorption  cell  also 
containing  only  H^O  and  N^,  The  absorption  path  is  uniform  with 
L - 100m,  T ’*■  300°K,  p = 0,  07  atm,  c = 0.0143  (case  9R),  and 
c = 0,  143  (case  IIMR). 

The  control  card  and  auxiliary  data  deck  input  for  four  radiance 
calculation  runs  is  shown  in  Figure  17  and  discussed  below  (numbers 
refer  to  indicated  card  groups  in  Figure  17). 

1.  The  calculation  is  performed  in  the  coupled  mode  with  the 
Lorentz  line  shape. 

2.  Since  the  absorption  path  is  uniform,  NAPATH  is  set  to  1 
and  MODEATM  is  set  to  UNIFORM.  Tlie  emission  path 
segment  is  also  uniform,  and  computation  of  the  unattenuated 
source  radiance  is  niade  wilJi  an  effective  value  NSPATM  1 
because  MODESRC  is  indicated  as  UNIFORM.  When  the 
coupled  path  calculations  are  made,  the  emission  path  seg- 
ment is  divided  into  10  integration  intervals  and  the  CG 
approximation  is  employt'd.  The  value  NSPATH  = 10  was 
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Figure  17.  Control  Card  and  Auxiliary  Data  Input  Deck 
for  ERIM  Hot-through- Cold  H2O  Spectra 
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determined  empirically  by  repeated  runs  with  increasing 
NSPATH  at  the  single  spectral  position  v = 3750  cm  ^ and 
noting  the  rate  of  convergence  of  the  calculated  radiance 
to  a constant  value. 

3.  The  spectral  region  from  3000  to  4400  cm  ^ is  covered 
with  56  equal  size  steps  (25  cm  ^ wide). 

4.  Only  H^O  is  considered  as  an  active  species  and  is  identi- 
fied by  the  name  H20  and  species  number  1.  The  combined 
band  model  parameters  for  H^O  are  used. 

5.  The  source  is  specified  as  a single  line  of  sight  with  the 
pTc  data  to  be  read  from  a card  deck.  Concentration  data 
for  H^O  will  be  taken  from  the  first  species  column  of  the 
input  pTc  source  deck. 

6.  The  pTc  data  deck  for  the  source  immediately  follows  the 
SOURCE  control  card.  Since  the  source  is  uniform,  only 
two  data  cards  are  required  to  specify  the  path. 

7.  The  absorption  path  is  identified  as  a single  line  of  sight 
with  the  pTc  data  to  be  read  from  the  following  input  deck. 
Concentration  data  is  to  be  taken  from  species  column  1 
for  this  run. 

8.  Compute  radiance  spectra  for  case  9R  with  the  CG  approxi- 
mation. 

9.  The  same  case  is  rerun,  but  the  attenuated  source  radiance 
is  computed  with  the  DR  approximation. 

10.  Case  1 IMR  is  run  by  redefining  the  species  column  in  the 
ab  sorption  pTc  deck  where  the  H^O  concentration  is  to  be 
obtained.  The  DR  approximation  from  the  previous  run  is 
retained. 

11.  Case  IIMR  is  rerun  with  the  CG  approximation. 
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This  job  run  will  result  in  lour  output  listings.  I lit  listing  lor  case 
9R  with  the  DR  approximation  is  given  in  Figure  18.  Ilie  eoniputed 
radiance  spectra  for  the  four  cases  are  compared  witli  experimental 
spectra  in  Figure  19.  A discussion  of  these  comparisons  is  given  iri 
Ref.  (13). 

B.  JS7— Jet  Engine  Plume 


This  example  illustrates  the  use  c^f  ATLES  for  a dc'tailed  study 
of  the  infrared  signature  of  a ji-t  - ai  r r raft  plume  in  both  the  2.  7-  and 
4.3-j.m  spectral  rogictns  and  for  a i rc  raft  - to- air  c ra  ft  and  space-to- 
aircraft  detection  geometries.  The  the rmodynaniic  properties  of  the 
plume  were  obtained  witli  the  General  Electric  code  SCORPIO  II  and 
represent  the  plume  of  a J - 57  engine  operating  at  full  military  power 
at  an  altitude  of  32,000  ft  (9.75km).  The  temperature  and  concentra- 
tion profiles  of  H^O  and  CO^  for  the  plume  are  shown  in  Figures  20 
and  21,  respectively,  and  are  tabulated  in  the  input  data  listing  of 
Figure  22.  The  plume  pressure  is  assumed  constant  tliroughoul  the 
plume  and  equal  to  the  ambient  pressure  at  9.75km.  A control  volume 
with  nozzle  radius  Rq  - 0.  35m,  length  L - 51m,  and  maximum  radius 

R - 1.  50m  is  used  to  contain  the  plume, 
m 

Botli  ll^O  and  CO^  are  considered  as  actixe  )jlume  and  atmos- 
pheric species  in  the  2.7-u.m  region.  Only  CO;,  is  considered  activi 
in  the  4.3-um  region,  '['he  calculations  are  carried  througli  with  a 
spectral  unit  of  cm  The  spectral  limits  m the  2 . 7 - am  region  are 

3000  to  4200  t in  ^ with  a step  size  ,'  \ - 25  cm  * (NSPFC  48).  In  the 

4.  3-um  region,  the  limits  are  2150  to  2400  cm  ^ with  \ 5 ^in  ^ 

(NSPEC  - 50).  Combined  banci  model  pa  rameters  are  used  and  selected 
for  a 25  t ni  spectral  resolution  in  the  2 . 7 - urn  reeion  and  5 tin  ^ 


7 - uin  re  u i i in  a nd  5 tin 


1'  . S.  Simmons,  I he  Aerosjxiee  (.or|).  , FI  Seeuntio,  l alil.  , private 
I oinnuini  i at  ion  (March  l'*75). 
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Figure  18.  Output  Listing  for  ERLM  Case 
9R  with  DR  Approximation 
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Figure  19.  Comparison  of  Computed  and  Experimental  Radiance 
Spectra  for  ERIM  Hot-th ru- Cold  H2O 

Measurements,  (a)  ERIM  9R;  (upper) 

experimental  source  radiance; (lower) 

experimental  attenuated  radiance: 

(upper)  computed  source  radiance; 

(lower)  computed  attenuated  radiance  (CG  approxi- 
mation);   compu'-ed  attenuated  radiance 

(DR  approximation),  (b)  ERIM  1 1 MR:  

experimental  attenuated  radiance; 

computed  attenuated  radiance  (CG  approximation); 

computed  attenuated  radiance  (DR 

approximation) . 
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resolution  in  the  4.3-jijri  region.  Besides  spectral  results,  in-band 
radiance  values  art'  calculated  for  two  bandpasses  in  each  of  the  two 
spectral  regions.  The  bandpass  limits  encompass  the  red-  and  blue- 
spikes  of  the  two  spectral  regions,  and  are  tabulated  in  the  input  data 
listing  of  Figure  22. 

Two  viewing  geonaetries  are  considered.  The  airc  raft-to-airc  raft 
geometry  assumes  a sensor  located  at  the  same  altitude  as  the  target 
plume  and  at  a horizontal  range  of  2km.  The  atmospheric  slant  patJi  is 
considered  as  a uniform  absorption  path.  The  s pac:e  - to- a ire  raft  geo- 
metry assumes  a satellite-based  sensor  70^^  off  the  target  zenith.  The 
atmosphere  is  assumed  to  terminate  at  100km  for  this  case.  The  por- 
tion of  the  atniosphere  between  the  plume  altitude  (9.  75km)  and  100km 
is  divided  into  10  intervals  of  approximately  each  mass  content.  Non- 
uniformities along  this  slant  path  are  handled  with  the  C'G  approximation. 

For  both  viewing  geometries,  the  atmospheric  H,0  c^oncentration  of 

(2 ) 

the  McClatchey  et  al.  tropical  model  is  used,  and  CO^  is  assumed 

uniformly  mixed  with  mixing  ratio  0.  000330.  A broaclsidc  (o  = 90°)  r 

aspect  is  assumed  for  both  cases. 

By  inspection  of  the  plume  profiles  and  the  magnitude  of  the  spec-  4 

tral  absorption  coefficients  for  H^O  and  GO^,  it  can  be  deduced  Lliat  the  f 

plurru;  is  relatively  optically  thin  throughout  the  spectral  regions  cem- 
sidered.  Tlie  integration  grid  for  the  emission  po'*tion  of  the  line  of 
sight  is  chosen  primarily  to  accommodate  the  degree  of  teiTipe rature 
and  concentration  gradient  within  the  plume.  The  calculations  are 
performed  with  an  integration  step  size  of  As  ^ 0.  1 Ojn  regardless  of 
where  the  line  of  sigh;  inti-rsecis  tlie  plunie.  Account  of  nonunilormi- 
ties  within  the  plume  and  that  result  from  tlie  coupling  of  tlie  plume  to 
the  atmosphere  is  made  with  the  DR  appro.ximation. 

Betore  complete  radiance  calculations  were  made,  an  extensi\e 
mapping  of  the  unattenuated  plume  radiance  at  the  single  spectral  posi- 
tion V - 3700  cm  ^ was  made  in  order  to  defini'  the  e.xtent  and  resolution 
of  the  area  integration  grid  that  would  be  required  in  c'rder  to  perform 
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accurate  Integrations  over  the  projected  area  of  the-  plume.  The  into-  | 

gration  grid  tabulated  in  the  inp>ut  data  listing  of  Figure  2.1  and  pictured  I 

in  Figure  23  was  finally  selected  and  verified  to  yield  integration  results  j 

accurate  to  approximately  two  percent  (at  3700  cm  ^).  | 

The  input  data  deck  for  the  final  runs  is  listed  in  Figure  22.  The  j 

first  run  covers  the  2.  7-p.m  spectral  region  and  the  horizontal  viewing  1 

geometry  case.  The  second  run  redefines  the  viewing  geometry  only.  j 

The  third  run  retains  the  space-sensor  geometry  of  the  second  run,  but  ^ 

redefines  the  spectral  region  and  active  species.  The  fourth  run  retains  i 

the  4.  3-gm  spectral  conditions,  but  reverts  to  the  horizontal  geometry. 

In  all  cases,  a spatial  map  is  made  for  the  two  bandpasses  considered. 

The  output  listing  of  Figure  24  is  limited  to  the  aircraft- sensor  geometry 
for  the  4.3-gm  region  and  blue-spike  bandpass  2375  - 2400  cm  The 
radiance  and  transmittance  spectra  for  all  of  the  runs  are  shown  in 
Figures  25  through  27.  Axial  station  and  isoradiance  contour  plots 
derived  (by  hand)  from  the  spatial  distribution  map  of  Figure  24  are 
shown  in  Figures  28  and  29,  respectively. 
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Figure  24.  Output  Listing  for  Jet-Plume 
Calculation  in  the  4.3-pm 
Region  with  Aircraft-Sensor 
Geometry  (sheet  2 of  2) 
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Figure  25.  Radiance  Spectra  in  2.7-|j.m 
Region,  (a)  Unattenuated 
and  Attenuated  Radiance  for 
both  Viewing  Geometries, 
(b)  Attenuated  Radiance 
(expanded  scale)  for 
Aircraft-Sensor,  (c) 
Attenuated  Radiance 
(expanded  scale)  for  Space- 
Sensor. 


TRANSMITTANCE  TRANSMITTANCE 


Figure  Z7, 


Atmospheric  and  Effective  Atmospheric 
Transmittance  Spectra,  (a)  Horizontal- 
Path,  2.7-pm  Region,  (b)  Slant-Path  to 
Space,  2.7-(j,m  Region,  (c)  Horizontal- 
Path,  4.  3-pm  Region.  (^)  Slant-Path  to 
Space,  4.3-p.m  Region.  T is  the  atmo- 
spheric transmittance  for  continuum 


radiation.  Xg  is  the  effective  atmo- 


spheric transmittance  defined  as  the 
ratio  of  attenuated  to  unattenuated  plume 
radiance . 
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VI,  CODE  DISTRIBUTION 


The  code,  band  model  parameters,  and  standard  atmospheres 
are  available  on  magnetic  tape  to  interested  users.  (A  blank  tape  must 
be  supplied  by  the  requestor.)  The  tape  consists  of  80  column  card 
images  blocked  at  60  cards  per  block,  and  can  be  supplied  either ;as  a 
9-track  tape  with  EBCDIC  coded  characters  or  as  a 7-track  tapejwith 
BCD  coded  characters.  Any  of  the  commonly  available  bit  densities 
can  be  accommodated.  The  tape  is  divided  into  17  partitions  separated 
by  end-of-file  marks.  The  first  partition  is  the  FORTRAN  source  deck 
for  the  code.'  The  next  10  partitions  are  the  band  model  parameter  sets 
of  Table  2.  The  last  six  partitions  of  the  tape  are  the  standard  model 
atmospheres  of  Table  1.  Each  of  the  10  band  model  parameter  parti- 
tions corresponds  to  the  card  deck  structure  of  Figure  12,  and  each  of 
the  six  atmosphere  partifions  corresponds  to  the  card  deck  structure 
of  Figure  13  (with  lUNIT  1),  An  input  deck  of  control  cards  and 
auxiliary  data  decks  corresponding  to  a limited  version  of  example  B 
of  Section  V will  be  supplied,  as  well  as  the  results  of  a job  run  made 
directly  with  this  deck  and  the  user's  tape. 

The  code  is  written  for  the  CDC  7600  computer  with  SCOPE2 
operating  system.  The  following  information  is  given  to  facilitate  the 
adaptation  of  the  code  to  other  systems. 

1.  Small  Core  Memory  Requirement:  125642g  words  (This 
includes  the  requirements  for  subroutine  SPLOT  and  the 
Aerospace  in-house  plotting  routines,  6307g  words). 

2.  Large  Core  Memory  Requirements:  203720g  words. 

3.  Computer  memory  must  be  preset  to  zero. 

4.  All  Hollerith  variables  and  constants  are  10  characters 
long. 

i 

5.  Input  and  output  devices  must  be  defined  to  be  consistent 
with  the  following  FORTRAN  statements: 


READ  (5,  format  number)  ....  card  input  (e.g.  , control 
card  deck). 

WRITE  (6,  format  number)  ....  printer  listing 

READ  (2,  format  number)  ....  band  model  parameters 

READ  (3,  format  number)  ....  standard  atmospheres 

Although  the  code  has  been  used  and  checked  on  many  test  cases, 
the  author  cannot  guarantee  that  the  coding  is  completely  error  free, 
and  in  fact,  would  appreciate  notification  of  errors  found  by  users. 

I • 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experinriental  and  theorotical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans* 
fer,  reentry  physics,  chernlcal  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high^power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos* 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  La^ratory;  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plaama  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging:  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal 
matrix  composites  and  new  forme  ot  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapone  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue -induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory;  Atmospheric  and  ionospheric  physics,  radia- 
tion frmn  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physice.  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionoepherc,  and 
magnetosphere:  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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